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ABSTRACT 


This woik amis to incorporate load chaiacteiistics in power flow and stability studies 
A computer code has been developed for modelling various kinds of loads It gives real 
and reactive powers of different loads for specified voltage and frequency This code 
has been linked with a power flow program the composite program with some addition 
facilities is named ’ LOADMOD 

LOADMOD consists of two main parts namely Load Flow Part (LFP) and Load 
Modelling Part (LMP) In LMP facility is given to the user to define a number of models, 
then parameters choiv^e of buses for load modelling to be applied, and distribution of 
load among defined models Compensation at the load buses is also modelled Real 
and reactive power loss in the feeders have been modelled, for which the initial real 
md the reactive power losses at the reference voltage and frequency are needed LFP 
interacts with LMP, i e LFP gives voltage magnitudes at different buses to LMP, LMP 
calculates P and Q and gives them as input to LFP and so on till voltages converge 

LOADMOD has been tested on NREB (357 bus) and IEEE (30 bus) systems with 
some assumed load mix at some chosen buses Substantial changes in the final voltages 
was observed with respect to the results for direct specification of real and reactive 
powers at load buses Facility has been provided in LOADMOD to iteratively inciease 
the loading by a fixed factor at specified or all load buses m the system, at the end of 
each LFP+LMP run tiU LFP fails to converge This failure to converge indicates that 
the system has reached the static stability limit Changes in the voltages have been 
studied with increasing loading for the IEEE 30 bus system in order to capture the 
buses prone to voltage instability 

Facility has been provided in LOADMOD for load and generator outages when the 
loading on the system is increased The results for such contingencies m IEEE 30 bus 
system have been studied 


VI 


Chapter 1 


Introduction 


1 1 Requirement of load modelling 

Electric power systems consist of generation transmission and load systems The per 
foimance of power systems is generally studied by simulations wherein mathematical 
equations of different parts are formulated, hnked properly and then solved These 
mathematical equations form the models for the corresponding parts of the power sys 
tern Till 1980 s generation and transmission systems were represented by mathemati 
cal equations of different degrees of complexity but loads were generally represented by 
constant real and reactive powers This representation was satisfactory in most cases 
of load flow and transient and dynamic stability analysis When the power systems 
became large and complex and the demand situation forced operation of the power 
network under heavily loaded condition, the simulated and measured values differed 
significantly especially with regard to voltage levels, and in some cases system collapses 
could not be explained as the simulation showed the safer side whereas the actual sys 
tern exhibited abnormal behaviour In such cases it becomes important to incorporate 
load behaviour in static as well as dynamic simulations Proper load modelling is thus 
important, especially in mid term and long term stability studies 
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1 2 Backgiound of this work 

The effect of load chai \cteiistics on dynamic perfoimance of powei system was known 
long back in 1930 s [1] riieii in 1960 s tlieie was work on incoipoiating voltage de 
pendence of loads and by 1982 much work was repoited on load representation [1] A 
problem with load modelling is, the loads are diverse and no single standard model can 
adequately represent different kinds of loads Furthermore, the actual load mix in an 
uer and load nixgmtudcs aie not easy to ascertain In mid 1980’s Electrical Power 
Research Institute of USA (EPRI) took a step to piomote woik on load modelling 
Ccncial Electric Company of USA was awarded a project for developing software for 
load modelling GE used the load models that weie available till date and developed a 
coniimic i ]) ickigc ( died LOADSYN [2] In the 1990’s IEEE task foico on load tepio 
scntation gave suggestions on standardizing some of the load models [4 5] There have 
been exlcnsivc discussions regarding various aspects of load modelling m the htratuic, 
[9 10, 11, 12, 13, 14 15, 16 17 18] 

1 3 Pxesent work 

This work is on the same lines as done by GE A package named "LOADMOD” has 
been developed in which the various models available have been incorporated and linked 
with a power flow program for the purpose of power flow analysis Facility has also been 
provided foi sialic slabdily siiidy 

Static load model parameters of different loads have been taken from [3] and from 
suggestions in [1, 4, 5] Computer code has been written for dynamic nrodcl of induction 
motor Classical model of induction motor was considered [6] 

Some load modelling concepts are discussed m Chapter 2 In Chapter 3 program 
development, the features provided to the user to study voltage profiles with increasing 
load and other features incorporated are discussed Test cases are considered and the 
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coiicsponding lesults -xie given in Chaptei 4 Cliaptei-5 concludes with discussion, 
sticngtlis and liinitations of this piogiam and suggestions for future woik 

Appendix A gives the computei code in FORTRAN foi the load modelling pail 
(LMP) of LOADMOD In Appendix B the mathematical backgiound of dynamic model 
of induction motoi is given, and its code is given in Appendix C The input and output 
files of LOADMOD aie given in Appendix D 


Chapter 2 


Load modelling concepts 


The concepts of load modelling and the ways and means of modelling are discussed in 
several papers and books [2 3, 4 5] In this chapter some of them and those applied 
to the present work are discussed 

Load modelling is considered in two mam approaches, namely Measurement based 
approach and Component based approach 


2 1 Measurement-based approach 

In measurement based approach measuiing devices are placed at the load buses that 
are to be represented These devices measure the variation in load active and reactive 
powers and the corresponding changes in voltage and frequency, (either due to inten 
tionally created disturbance or naturally occurring) From these measurements, the 
aggregated parameters of the model for the total load at the bus are derived 

The advantage of this method is that it gives the aggregated load model parameter 
of the actual system directly But it has the following disadvantages 
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« Cost of the measuring equipment may be significant 

• Tlieie IS need to make measurements at all buses which aie to be represented, 
because the loads at different buses generally have diffeient load patterns and 
mix If two buses have similar loads then measurements at one bus are sufficient 
and the same parameters can be apphed to other 

• As the load mix changes at different times of day and seasons continuous mea 
siuement under these varying conditions is required 

2 2 Component-based approach 

In component based approach detailed modelling is done for all or important compo 
nehts of loads Then either the model parameters are aggregated to get one load model 
01 they are used independently In this work they are used independently i e , each 
mathematical model for a component at a bus is activated by giving voltage and fre 
ciuency inputs which give the real (P) and reactive (Q) power outputs Then the P and 
Q outputs of all models aie added to get the total load at the bus 

This method overcomes some of the disadvantages of the previous methods Thus, 

• By changing the mix of the load, the load characteristics may be changed to suit 
the particular load mix at a given time or to suit seasonal changes 

• No hardware is required, only software is needed to be implemented 
The disadvantages of this approach are, 

• As the loads are diverse, data collection regarding the number of load components 
and load composition is the main problem 

• New load models are to be worked out and implemented when new or redesigned 
equipment is installed 
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• Moie computational time is needed 

Foi Indian conditions where it becomes difficult to make measurements on the sys 
tern measurement based approach is difficult to be tested and implemented The 
component based approach is better suited With sufficient statistical data on how 
the load mix changes in different areas at different times and seasons this approach is 
expected to give satisfactory results with less cost and problems aa compared to the 
measurement based approach 

Foi the survey of load mix and their changes with seasons EPRI used consumer 
bill data Similar studies can be done to obtain this data and implement it for Indian 
conditions 

The load model is of two types, namely, static and dynamic 

2 2 1 Static load modelling 

In the static load model the real and reactive powers at any instant of time are expressed 
as functions of bus voltage and frequency at that instant 

All load components can be represented by static models Examples of load com 
ponents are lighting, heating, motor loads , air conditioning, arc furnace electrolysis, 
etc 


Static models are of many types [2 3, 4, 6] The most common among them are 

• Constant impedance load model 

• Constant current load model 

• Constant power load model 
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Polynomial load model 

In polynomial load model load active and reactive poweis are expressed as 

= +B(^) + C1 (2 1) 

Q = QAO(^? + E(^) + F] (2 2) 

Vo Vo 

Similarly, P and Q can be expressed as functions of frequency with (///o) replacing 
{V /Vq) in the above equations and different constants 

Heie Po and Qq are the values of P and Q at nominal values of V and / i e at Vq 
and /o P and Q aie the alues of active and reactive power of the load at any V and f 

This model is also called the ZIP model [6] as the earlier three models can be achieved 
from this by setting some of its constants equal to zero 

For example, in equations 2 1 and 2 2, if B, C, E and F are zero then it becomes 
a constant impedance model Similarly other models can also be achieved Hence 
polynomial model is a combination of the earher three models and can give load P and 
Q in between the results of the constant Z, I or P models 


Exponential load model 

In exponential load modeUing, load active and reactive powers are expressed as 

p = Poiy/ (2 3) 

Q = «o(h)0 (2 4) 

Similarly, they can be expressed for frequency dependency Sometimes two or more 
terms with different exponents are used 

If some equipment behaves in a particular way for a particular range of voltage and 
in some other way beyond that range it can be easily incorporated in the exponential 
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model scheme by having two exponents one for each iduge as 

P=^Po{^f^forVi<V<V2 

yo 

P = Po{^f^forV2<V<V, 
yo 

The numbei of ranges can be more than two This feature can be incorporated for the 
polynomial model also but it would require more number of model parameters 

In this work the exponential model is used for static modelling The data for various 
models were taken from Taylor [3] 

In the equations beiow paiametei nm represents the fraction of the total equipment 
load that is motor Equation 2 5 and 2 6 below are used when the equipment is either 
fully motor or not at all whereas equation 2 7 and 2 8 are used when the equipment 


contains some part as motor as in heat pumps etc 
For nm=0 0 or nm=l 0 

P = Pol(h)'- {Ifi] (2 6) 

Q = Ool(h)« (j-fn (2 6) 

For 0 0 < nm <10 

P = nmP„l(h)f (l)P,\ + (1 - nm)P,[(~Y‘~ {{)’’■ »] (2 7) 

ko Jo *^0 Jo 

« = «m<2„[(h)« (■f)«'l + (l-nm)Qo!(h)e "(k)0/ -] (2 8) 

VQ Jo yo Jo 


In the above equations, 

P = Bus real power as given by the model 
Q — Bus reactive power as given by the model 
Po = Initial real power at rated voltage and frequency 
Qo — Initial reactive power at rated voltage and frequency 
= Voltage dependent exponent of real power 
Qy = Voltage dependent exponent of reactive power 
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Pf = Fiequency dependent exponent of real power 

Qf = Fiequency dependent exponent of reactive power and 

Pv m Qvnm Pf m and Qfnm are the voltage and frequency dependent exponents of 
active and leactue povers respectively of non motor part of load 

The piocedure followed in this work is that the load at the bus is divided among 
the defined models (foi which the parameters of exponential models are available) i e 
pQ foi each model is given From Pq and power factor reactive power Qq is calculated 
The model is activated by the V and / inputs where V input is obtained from Load 
Flow Piogram and / is calculated fiom the frequency regulation characteristics of the 
system Equations 2 5 2 6 2 7 and 2 8 are used to calculate P and Q from the input 
V and / values 

Similarly, this procedure is repeated for all the models and their corresponding real 
and reactive powers are obtained Individual P’s and Q s are then added to get the 
total P and Q of the load at the bus for the given V and / Then the feeder real and 
reactive power losses are added to P and Q which completes the P and Q calculations 
at the bus for given V and / Because of varied voltages the line P and Q losses also 
vary The implementation of this is discussed in the later part of this chapter 

In EPRI LOADSYN, the individual load model parameters are aggregated and one 
load model is obtained m which the compensation at the bus is also mcluded In this 
case when the aggregated load reactive power is equal to the compensation the Qq in 
the equation 2 6 and 2 8 become zero This causes a problem since, the bus reactive 
power is calculated to be zero To avoid this problem LOADSYN normahses the reactive 
power formula to Pq instead of Qq In the present work as each load model is treated 
independently and at the end compensation, if any is added, the above problem does 
not arise 

With good models and accurate load mix representation, static load models, when 
used for static or dynamic simulations give better results than the the conventional 
constant P, Q models When LOADSYN was tested on Ontario hydro system with 
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onl}- static models significantly improved results weie obtained [2] 

When static models are used for large dynamic loads such as induction motors their 
dynamic behaviour may not be captured accurately The static models of such loads 
gi.'e inaccurate results It then becomes necessary to accuiately model the dynamic 
behaviour of the equipments 


2 2 2 Dynamic load modelling 

In the dynamic load model, the real and reactive power at any instant of time are often 
expressed as functions of bus voltage and frequency at the previous instance of time m 
the discrete time representation of system dynamics Alternatively in the continuous 
time formulation an appropriate differential equation may be used for load dynamics 
The components in the power system that require dynamic modelling are 

• Induction motors 

• Thermostat control of space heating and refrigeration equipment 

• Load tap changing effect of transformers 

• Over current protection of big induction motors 

• Synchronous motors 

• Arc extinction and restart of discharge lamps 

The extent to which these need to be modelled depends on the type of study to be 
done Generally, since most dynamic loads are induction motors, their modelling has 
been discussed extensively [2 5 6] 

In the present work, classical induction motor equivalent circuit is taken and the 
swing equation (differential equation) is solved using Runge Kutta method The math 
ematics of this is given in Appendix B and its computer code is given in Appendix C 
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If all the dynamically behaving equipments m the system are modelled propeily 
then accuiate lesults will be obtained and probably all cases of voltage collapses would 
be explained properly 


2 3 Modelling compensation 

C eneially at load buses reactive power compensating equipment is present This com 
pensation vanes with change in voltage at the bus 

AT 1 0 p u voltage i e for Vi = 1 Op tt at the bus the reactive power injection Qi 
at the bus for fixed capacitive shunt compensation is 

Qi = Vi^y (2 9) 

When the bus voltage changes to V 2 we have 

<52 = V2^y (2 10 ) 

From equation 2 9 and 2 10 

Q2_V^ 

Qi U' 

< 3 . = 

In equation 2 9 Vi = lOpu so that Qi = y 

Qi = Vj’Qi (2 11) 

Equation 2 11 is implemented to get the compensation at any voltage 


2 4 Modelling real and reactive feeder losses 


Let the initial real power loss in the feeders at 1 0 p u voltage at the terminal of the 
load be specified 
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Let R = resistance of distribution line This is a fixed quantity Then the real power 
loss on the distribution line is 


Li = 


(2 12 ) 


01 




(2 13) 


Wheie Pi and Qi are load real and reactive powers at 1 0 p u voltage at the terminals 
of the load 

Let the losses at the changed voltage V 2 be L 2 
Then 


52^ „ Pz' + Q2^ 


L2 = ^R = 




R 


■'2 »'2 

Pa and Q 2 are obtained from static load modelling Prom equations 2 13 and 2 14 

Li 

L 


(2 14) 


u 


Li .Pa^ + Q2' 


•( 


) 


Va'^Pi' + gi '*- 
Equation 2 15 is implemented for the losses at any voltage 


(2 15) 


Similarly, if the imtial reactive power loss at 1 0 p u voltage at the terminal of the 
loads or the reactance of the distribution hne is known, Qio,, can be calculated at 
any voltage and the same equation 2 15 can be applied, where, m this case L 2 would 
represent Qi ss under various conditions of load This modeUing of real and reactive 
power loss by equation 2 15 has been implemented in LMP Flexibihty is given to the 
user to choose an option to use this facility 
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In the system because of the diops m R and X of the distribution lines the actual 
volt !{,( s it the lo rd component is less than that at the bus If Pi and Qi aie the ical and 
le lotive poweis of load it 1 0 p u voltage it implies that the voltage at the component 
is 1 0 p 11 The bus voltage will not be 1 0 p u because of line chops Howcvci, in 
calculating P\ and Qi the chstiibution line diops have been neglected Hence the actual 
P and Q of loids aie not cciual to Pi incl Qi calculated neglecting line chops An 
itciative piocccluie can be developed to accuiately incorpoiate line drops and losses 
This, howcvei has not been done since, if the links aie strong enough, reasonably 
accuiate lesults aie expected without lesoiting to iterations However m future this 
can be taken up to modify LMP to get more accurate results The computational 
complexities will, however increase 

Iiidii has luge igiicultural loads A major component of those loads arc incUiction 
motors, fed from long distribution lines, because of which line drops are high This 
implies that the appioxiniate foiinulae derived earlier for line loss calculation may be 
unsatisfactoiy foi agiicultuial loads Hence accuiate modelling of load as well as chops 
m lines may be required m the Indian context 

In (he pioscnl woik, s(a(ic model paiamctcis aie taken from Taylor [3] Tliough 
these pu imetcis aie lelcvant loi Indian conditions, there aic aspects of Indian loads 
which are rmiciue Poi instance, the large scale use of small uninteiinpted power supply 
units by comineicial and domestic users, stems from fiequent power breakdowns and 
annormced power cuts Load modelling should incorporate such units There may be 
ispccts m dynamic loads ;which aio pcculiai to India There is Unis need for considoimg 
the charactei of loads in India in detail and developing new and lelevant models both 
static and dynamic The scope of the present thesis however, does not include such a 
study 

Load modelling when observed horn these angles is a vast field and gives much scope 
for further work This component based appioach would be most successful when load 
statistics aio properly known It is cost effective By simple vaiiation m input data files 


Hk computci piogi uii foi load modcllmg can be made suitable to any system 

In the next chaptei the flexibilities piovided in LOADMOD and the imiilementation 
of the piogi un aie discussed The mam piogiam with all the data files lequucd is given 
in a floppy disc with this woik 


Chapter 3 

Programming details 


3 1 Introduction 

lu this chapiei the load modelling piogiam and the details of its linking with the load 
flow piogn .111 me discussed The code wiitten in FORTRAN 77 is given in Appendix A 

The load flow pait has one input data file ’Ifdata’ whcieas load modelling pait 
has scvcial input data hies whose details aie explained in the following sections The 
LOADMOD with all input data files and the code aie given in a floppy disc with 
this woilv In the file named ’iiifoiination’ all the data files lequiied aie listed An 
infoiination file is given with each input data file which has the same name as the data 
file with an extension ’ inf Foi example, the data file ’gen outage has infoimation file 
named gen_outagc mf ’ 

Ihe load modelling pait is divided into foui modules whicn aie appiopiiately linked 
with the load flow pait The LFP used is a Newton Raphson load flow piogiam in polai 
CO oidinitcs using spusity Lccliinquo 

The details of the modules and their linking with load flow progiam aie discussed 
in the following sections The flow chut of the LOADMOD is shown in Figuic 3 1 


1/1 
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3 2 Module - 1 

Module 1 essentially serves two basic purposes First all the variables used in LMP are 
declared Secondly some of the input and output files required foi LMP are opened 
This module appears along with the declarations of the Load Flow Part The details of 
some of the main variables which are used for linking of LFP and LMP are discussed 
m detail in the next modules 

3 3 Module - 2 

This module stores the compensation, if any, at the load bus in an array ’qq(bus mim 
bei) which is subsequently made use of in the Load Modellmg Part This module is 
placed in the Load Flow Part where the data for the load buses is read 


3 4 Module - 3 

This module is placed after all the data reqmred for the Load Flow Part is read This 
module reads the data needed to provide flexibility to LOADMOD The input files 
which are read by this module are as follows 

3 4 1 File ’itrJimit’ 

This file sets limits in the iterations of 

• LFP 


• LFP+LMP convergence 
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Flexibility has been piovided to either keep the value of voltage tolerance constant 
niespective of numbei of iterations oi to change the toleiance after some required 
numbei of iterations of LFP+LMP runs 

3 4 2 File ’gen_outage’ 

This file stores the information regarding generator outages which are as follows 

• Whether generator outage is required or not 

• Number of generators 

• Generator bus numbers 

• The load increment iteration after which the generator is to be taken out 

3 4 3 File ’load_outages’ 

This file contains information similar to that in ’gen-Outage’ 

3 4 4 File ’out_plot’ and ’out_plot_names’ 

These files facilitate LOADMOD to generate individual output files for each bus in 
which the bus voltage and P and Q values at the bus throughout the different states for 
the system (i e , system states at the end of each LFP+LMP convergence) are written 
for easy checks and plotting The choices to be given in these files are 

• Option for using this facility 

• Selection of all or a hmited number of buses If all option is chosen, then LOAD 
MOD generates output files for each bus with a default file name ’o bus number’ 
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for which it uses file outJile-Qames to store the default generated file names If 
limited option is chosen the bus numbers and the required file names are to be 
specified in file out_plot Data foi all the options are read from this file 

3 4 5 File ’load_mix_change’ 

This file provides the option for changing the load mix from one state of the system to 
the next Once this facility is opted for then the following data are read from the file 

• buses at which the load mix change is required 

• the models where the loads change 

• the amount of load mix change 

3 4 6 File ’react Joss’ 

The loads at the buses are fed through distribution feeders because of which there is 
a reactive loss in the fines This Q loss depends on the square of the current which 
in turn depends on the apparent power of the load The loads change with change in 
bus voltages For modelling Q loss and taking it into account while caBculating the 
Q injection at a bus the data required is the reactive power loss on the distribution 
feeders assuming 1 0 p u voltage at the terimnal of the load This data i e the bus 
numbers and the Q loss at the bus are read from this input file 

3 4 7 Files ’load_modelling’ and ’load_composition’ 

For modelling real and reactive losses in the feeders at the bus equation 2 J5 is used In 
this equation SI represents the total apparent power of the loads at 1 0 p n voltage (at 
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the teiminals of the load) This calculation of Si at each bus is done here and stored 
in an airay si (bus number) for which the inputs requited are 

• Bus number at which load modelling is apphed 

• The total real power of load at these buses 

• The load composition at these buses 

The first two requirements are read from file load-modelling and the third from 
file load-Composition Then the load model subroutine model is activated which 
implements the load model characteristics and gives the total real power of the modelled 
load in variable pH’ and reactive power m variable ’q’ These values of pll and ’q 
are stored in a complex array ’si (bus number) The details of the subroutine model 
aie discussed in a later section 

At the end of this module two variables are used i6 and i2 At the end of each 
LFP+LMP run the execution is transfered to statement labeled 899 and variable 
’i2 keeps track of the iteration number i e , the number of iterations executed in the 
execution of the LFP+LMP run for one setting of load When the difference in bus 
voltages of present and previous LFP+LMP runs is less than the tolerance limit this 
process is stopped and the system power level is increased for next state of the system 
and the execution is transfered to statement labeled 8980’ This is called load increment 
iteration and is tracked by variable i6 


3 5 Module - 4 

This forms the major part of LMP m which 

• Data transfer between LFP and LMP occurs 

• Load characteristics are implemented 
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• Real and reacti've losses are modelled 

• Options specified in Module 3 are implemented 

• Sequence of operation of LFP and LMP are controlled 

This module is placed just at the end of the Load Flow Program The following steps 
describe the operation of this module 

3 5 1 Step - 1 

The buses at which loads are to be modelled are read into an array bn and then 
corresponding real powers into an array ’bpwr 

3 5 2 Step - 2 

If the option for load mix change is chosen then the following data is read from the file 
loadjmix-change’, for that state of the system 

• The buses at which the load mix is to be changed 

• Model numbers at which the loads are to be changed for each bus 

• The change in the model load for each bus 


3 5 3 Step - 3 

Each bus is considered one by one and the load model characteristics are implemented 
taking into account changes in the load mix and load outages as described in the fol 
lowing points 
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• The bus voltage is tiansfered from LFP (variable x’) to LMP (variable v ) Fre 
quency / and time step t which are read from file load-composition are kept 
constant In this work the frequency changes are not implemented as they are not 
as important as voltage changes This can be done in a straightforward manner 
if the frequency the frequency regulation characteristics of the system is imple 
mented When used with a stability program the values of / and t are fed from 
it 

• With the above information the subroutine model’ is activated which implements 
the load model characteristic and real power loss modelhng Subroutine model’ 
gives modelled real and reactive powers (including the real losses related to the 
bus) as its output m two variables p and ’q The details of the subroutine 
’model are dealt m detail later 

• The variables p’ and ’q’ are transfered to variables pload’ and qload’ which are 
used in LFP for real and reactive powers of the loads 

• The compensation at the load bus (variable ’qq* m module 2) is modelled using 
equation 2 11 and added to variable qload’ 

• Reactive power loss m the feeders is modelled using the equation 2 15, and is 
added to the variable ’qload’ With this the calculation of P and Q at the bus is 
complete 

• The LFP gives its real and reactive power injection outputs m the variables ’pinj 
and qinj’ The output of the modelled load at the specified buses is transfered 
back into ’pinj’ and ’qinj’ 

• The voltage values and real and reactive powers at the buses where the load 
modelling is apphed, are written mto a fiile ’itr out2’ for each LFP+LMP run 
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3 5 4 Step - 4 

This step checks for the convergence of LFP+LMP runs The actual operations m this 
step are as follows 

• The voltage of all the buses obtained from LFP are stored in an array xll 
Prom the second iteration onwards (of LFP+LMP run) the present voltages are 
compared with the previous iteration voltages ( xll ) and the maximum value of 
the diffeience is stored in variable xx3 

• The value in xx3 is compared with the voltage tolerance limit If it is more 
than or equal to the limit then the control is transfered to the statement labeled 
899 for repetition of another LFP+LMP iteration If the tolerance is satisfied 
LFP+LMP convergence is achieved and further options are implemented 

• When the system is over loaded (implying decrease in stability margin), or the 
toleiance limit is very small a facility is provided to change the voltage tolerance 
limit to some other value after some iterations of LFP+LMP runs 

3 5 5 Step - 5 

When LFP+LMP runs converge the bus voltages and the real and reactive power 
injections at all buses are written into the file ’itr outl’ At the end of LOADMOD, this 
file contains information regarding V, P and Q at all buses for aU states of the system 

3 5 6 Step - 6 

When the LFP+LMP conveiges for each state of the system, the voltage and real and 
reactive powers at the bus and total real and reactive powers of the system are written 
into the corresponding files for the buses as per the option in the data file out_plot’ 
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Thesp output files aie used to observe the voltages at the bus for different loadings of 
the system 


3 5 7 Step - 7 

Depending on the options in the data files gen outage and load-outage the generators 
and loads are taken out For load outages P and Q of the loads are made equal to 
zeio For generator outages PGEN PMAX and QMAX are made equal to zero and 
the variable IFLAG used by LFP to distinguish generator/load /slack bus is changed 
to make it a load bus Hence from the nex^' iteration this bus is treated as load bus 
with zero P and Q loads The convention used by the LFP for IFLAG’ is 

• IFLAG=1 for generator buses 

• IFLAG=2, for load buses 

• IFLAG=3, for slack buses 

• IFLAG=0, for no connection to the bus 


3 5 8 Step - 8 

For the first state of the system i e when ’i6 is equal to one, the following information 
IS read from the file ’load-increment 

• information regarding the options for other states 

• the factor of power level increment for the coming states 

• the buses where the powers are to be increased 
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3 5 9 Step - 9 

If the power inciement option is chosen then control is transferee! to the statement 
labeled 8989 i e, to the end of Module 3 The LFP+LMP convergence for the next 
state of the system is taken up 

Thus after some power increments of the system the voltages at buses start de 
cieasing and the number of iterations taken by load flow as well as the LFP+LMP to 
converge, increases This operation is continued till the load flow fails to converge i e 
the system reaches the static stability limit 

3 5 10 Subroutine - ’model’ 

The sequence of operations m this subroutine are discussed below 

• The bus numbers at which the load modelhng is to be applied is transfered to the 
subroutine ’model 

• From the input file ’comp.code’, the information regarding the number of models 
defined are read into a variable no_moder 

• The information regarding the load composition (real power at the bus divided 
among the defined models) is read from the input file ’load_composition’, into an 
array ml(model number)’ 

• Elements in the array ml’ are modified for 

— load increments 

— load outages 

— load mix changes 

• Real losses in the feeders catered to by the bus are read from the data file 
’load-comp osition’ 
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• Only the model chaiacteristics at that bus are implemented using a subroutine 
named modelchar , the details of which are dealt latei The information passed to 
the subioutine modelchar is 

— real power of each model (Pq of each model) 

— voltage at the bus 

— fxequency at the bus 

— time step (used when linked with stability program) 

The subroutine ’modelchar gives the modelled real and reactive powers of each model 
in the arrays realp and ’reactp All the elements m the arrays ’realp’ and reactp are 
added into the variables p’ and q respectively which give the total modelled real and 
reactive powers of all the models at the bus The real power losses in the feeders at the 
bus are modelled using equation 2 15 and are added to the variable p’ This completes 
the calculation of real power at the bus i e , p injection at the bus The same procedure 
IS repeated for all buses to get the corresponding modelled real powers 


3 5 11 Subroutine - ’modelchar’ 

In this subroutine, 

• Option for static (exponential or polynomial) or dynanuc modelling of each model 
IS read from file ’s_d’ 

• Parameters required for different exponential models are read horn file ’compxhar’ 

• The real power and power factor of each model is used to calculate the reactive 
power 

• Equations 2 5, 2 6, 2 7 and 2 8 are used to obtain the modelled real and reactive 
power loads of each model Facility is given to opt for a polynomial model, but 
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has not been de\ eloped in the current program In future work the parameters of 
the model aie to be fed in the file comp_charp’ and the mathematical equations 
aie required to be written 

Code IS developed for dynamic modelling of induction motor and is given in Appendix 
C This appears in the form of a subroutme dynmodel m LOADMOD Details of 
dynmodel are discussed below The real power input of the motor model may consist of 
motors of different range and each range may have a number of motors This information 
is lead in subroutme modelchar from the file d-d’ and supplied to the subroutine 
dynmodel’ The outputs from this subroutme are the modelled real and reactive powers 
of the motor When used with stability program the dynamic model of induction motor 
requires voltage, frequency and time step These should be supplied by the stability 
program 


3 5 12 Subroutine - ’dynmodel’ 

• The parameters of the motor of one range read from the input file ’d_d are 

— rated power 

— fraction of the rating at which the motor is used 

— power factor 

— Rs Xs, Xm, Rr and Xr of the equivalent circuit (Figure B 1) 

— inertia constant, H (MJ/MVA) 

— A, B and C values used in equation B 8 

• The value of a;* and Urn are calculated 

^base ~ 2 OlT fbase} fbase ^0 0 Bz 

— u = 2 Ott/s, /, = actual frequency 


— Um at full load = 013(1 0-slip) 
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The values of ujm at full load (oJmoi) are stored m variable ’wmOl 

• Full load torque is calculated 

- tmO = power/wmOl 

o The paiameters of the Thevenin equivalent circuit (Figure B 2) Re and Xe are 
calculated using equation B 2 

• Thevenin equivalent voltage is calculated using equation B 1 

• Runge Kutta method is used to solve the differential equation (equation B 9) and 
the value of utm is obtained 

• the value of tm for this value of Um is calculated using the equation B 8 appendix 
B 

• From the values of aim and tm the modelled real power is obtained The reactive 
power of the motor is calculated from this real power and the power factor 




This value of Um becomes the initial value for the next iteration of the model 



Chapter 4 


Results and discussion 


111 Lius clupLei the applicaLioii of LOADMOD to study the load flow and static stability 
of two systems is considered The first system is the system under the jurisdiction of 
Northern Regional Electricity Board (NREB) of India It is a large system with total 
load of about 17 000 MW The second system is the IEEE 30 bus system The aim of 
this study IS to do the load flow analysis with static load modelling and to capture the 
bus voltages with increasing system loading The details of the study are given below 


4 1 NREB system 

4 11 Load flow study for base load condition 

In the NREB system load modelling is applied at eleven buses of 220 and 132 K V level 
The details of the NREB system are shown in Table 4 1 [7] 

The buses chosen for load modelling and then real power loads at nominal values of 
voltages are given in Table 4 2 

The real power loads at these buses are divided among the defined models at each 
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Table 4 1 Details of NREB system 


Number of total buses 

375 

Transformers 

155 

Lines 

438 

Z loads 

80 

PQ loads 

166 

Generators 

70 

Slack bus 

l(Bus No 183) 


Table 4 2 Buses modelled m NREB system 


BUS NUMBER 

P_LOAD(MW) 

C 

o 

Q 

115 

121 08 

1 0 

134 

70 94 

10 

135 

19 87 

10 

136 

35 00 

10 

152 

24 59 

1 0 

158 

41 62 

10 

181 

202 48 

10 

243 

226 76 

10 

259 

247 77 

10 

267 

240 66 

10 

282 

317 09 

10 
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Table 4 3 Results at the end of one LFP+LMP iteration 


Bus No 

Volt(p u) 

Pdoad 

QJoad 

115 

8744 

110 6101 

25 4810 

134 

8382 

63 1387 

20 2399 

135 

8050 

17 2694 

3 1093 

136 

8472 

31 3541 

6 1463 

152 

8680 

22 3608 

9 7205 

158 

8582 

37 5732 

8 7956 

181 

9636 

197 1735 

85 6235 

243 

9712 

221 9774 

13 8056 

259 

9842 

244 9408 

22 4375 

267 

1 0033 

241 2382 

39 2232 

282 

9113 

297 7038 

28 2234 


bus and LOADMOD takes these inputs as the nominal power Pq, of each model at 1 0 
p u voltage At the end of one LFP run, the voltages at all buses are calculated and 
aie fed to LMP LMP calculates the values of real and reactive load powers, Pioad and 
Qi ad respectively at the buses specified where the load modelling is applied The bus 
voltages and the real and reactive power injections at the end of one LFP+LMP run 
are shown in Table 4 3 

From the above results it is observed that, at bus number 115, the Pioad calculated 
by LMP IS 110 61 MW instead of the specified value 121 08 MW (Table 4 2) and Qioad is 
as calculated by the load models Similarly, at the other buses shown above the P and 
Q loads are calculated using the voltages at the respective buses The Q load calculated 
mainly depends on the load mix data provided, i e , the division of the real power load 
among the defined models 

These real and reactive power loads are given fed to the LFP for another LFP+LMP 
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Table 4 4 Results at the end of LFP+LMP convergence 


Bus No 

Volt(p u) 

PJoad 

QJoad 

115 

9079 

113 3142 

38 1916 

134 

9225 

67 0958 

26 8317 

135 

9021 

18 5186 

6 1311 

136 

9117 

32 8464 

10 6632 

152 

9424 

23 5933 

11 0906 

158 

8851 

38 3063 

13 3028 

181 

9774 

199 1684 

89 0539 

243 

9846 

224 1940 

24 3254 

259 

1 0046 

248 5980 

23 4187 

267 

1 0178 

243 7789 

33 0406 

282 

9624 

308 7208 

61 2777 


run The voltage tolerance limit for convergence of LFP+LMP runs is taken to be 0 005 
p u, and is increased to 0 01 p u after 5 iterations Thus, when convergence is achieved, 
the difference in voltage magmtude at any bus for two consecutive iterations is less than 
the tolerance limit 

LFP+LMP convergence was achieved in 4 iterations and the end results are given 
m Table 4 4 

From these results, it is observed that when the bus voltages change the loads change 
(which is not taken into account in the conventional constant P, Q model), and as the 
loads change the bus voltages again change This is taken into account in LOADMOD 
LFP+LMP runs by exchanging the voltages and loads and they finally converge to the 
proper steady state values The amount of change in the voltages and loads depends on 
the load composition The final values of voltages and the real and reactive powers at 
the modelled buses depend on the load mix data and the load characteristics available 
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4 12 Static stability study by gradually increasing the system 
loading 

When the LFP+LMP conveiges and the operating point for that state of the system, is 
obtained the total system power level is increased by a factor of 1 01 i e the following 
quantities are multiplied by the factor 1 01 

• P 

gen 

• Pioad and Qioad at the load buses which are not modelled 

• Po of all the models at the buses which are modelled 

The same procedure as of earlier LFP+LMP runs is followed for getting the load flow 
for this case Thus the system power level is incremented systematically till the loadflow 
program fails to converge indicating that the system reached the static stabihty limit 
Voltages P and Q loads, at the buses modelled at a state prior to failure of LFP are 
shown in Table 4 5 

The voltage, Pioad and Qioad at the specified load buses (data file out_plot’) at 
different system loading are stored in separate files The voltages at bus 115 are shown 
in Table 4 6 

The variation of voltages at buses 115 134 and 135 with the total system real power 
IS shown in Figure 4 1 

To capture the voltages close to the stabihty limit, load increment was done with 
decreased increments around this point The resulting voltages at buses 115, 134, 135 
and 158 are shown m Figure 4 2 

By modelling the loads at all buses with accurate load mix data the voltages and 
loads around static stabihty limit can be captured accurately, which can be used to 
study the problematic buses 


I 
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Tiblo 4 5 Results of tlio state ])iioi to the failuic of LFP 


Bus No 

Volt(p u) 

PJoad 

QJoad 

115 

8406 

117 6417 

49 3039 

134 

8394 

68 8664 

28 8229 

135 

8118 

18 9144 

7 6729 

136 

8183 

33 4713 

13 6792 

152 

8648 

24 3075 

10 4655 

158 

7895 

38 9878 

15 4109 

181 

9623 

214 4887 

100 9493 

243 

9558 

238 9556 

96 9486 

259 

9571 

261 5258 

87 0234 

267 

9837 

258 9569 

33 2178 

282 

8907 

319 5829 

140 8106 


NREB SYSTEM 



Figiuc 4 1 Voltages at buses 115 134 and 135 with system load mciement of 1 01 
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BUS 115 


16878 410 4180 284 908 113 314 38 192 


42 462 


3 17210 053 4307 711 892 114 410 45 443 


4 17378 092 4370 181 884 115 012 47 395 


5 17547 734 4432 576 877 115 659 48 628 


4494 475 870 116 323 49 348 


7 17888 189 4583 042 862 116 902 49 994 


853 117 416 49 889 


9 18228 732 4745 934 \ 841 117 642 ' 49 304 


(This data is stored by LOADMOD in the file name 
specified by the user in file out_plot ) 


A System state number(first one is the base case) 
B Total P Joad of the system in MW 
C Total QToad of the system in MW 
D Voltage at the bus m p u 
E P Joad at the bus in MW 
F QJoad at the bus in MW 
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Total system real power load in MW 


Figure 4 2 Voltages at buses 115 134 135 and 158 with system load increment mod 
ificd to capluie the stability limit 

4 2 IEEE 30 Bus system 

4 2 1 Load flow study 

IFEE 30 bus system is a bench maik system generally used for powei flow and strbihty 
studies Theie aie othei bench mirk 4 bus system, 10 bus system, 75 bus system etc 
In this study 30 bus system is taken [7] Table 4 7 gives the base data about the numbei 
of buses, transformers, generators, etc of the IEEE 30 bus system 

Load modelling is applied at nine buses where the loads are high The details of the 
buses and their real power at nominal values of voltages are given m Table 4 8 

The leactive loss in the fcedeis at 1 0 p u voltage at the terminals of the load which 
foim tlio input foi leactivo loss modollmg aic assumed to be as shown in Table 4 9 

For the base case, LFP+LMP conveiged in four iterations The bus voltages and the 




Table 4 7 Details of IEEE 30 bus system 


Number of buses 

30 

Number of transformers 

4 

Number of lines 

37 

Number of load buses 

24 

Number of Generators 

6 

Slack bus number 

l(Bus No 1) 


Table 4 8 Buses modelled lu IEEE 30 bus system 


BUS NUMBER 

PXOAD(MW) 

Voltage(p u) 

8 

58 

10 

9 

112 

10 

11 

76 

10 

12 

22 8 

10 

14 

62 

10 

15 

82 

10 

17 

90 

10 

21 

17 5 

10 

30 

10 6 

10 
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Table 4 9 Assumed initial reactive losses at the modelled buses 


Bus no 

QJoss(MVAR) 

8 

0 1 

9 

0 35 

11 

0 15 

12 

1 2 

14 

0 2 

15 

03 

17 

07 

21 

03 

30 

05 


coirespondmg real and reactive powers at the modelled buses are shown in Table 4 10 
The loads at the modelled buses are as calculated by the LMP, in which the real and 
reactive losses and the compensation is also modelled 


4 2 2 Static stability study 

The load on the system was incremented by a factor of 1 1 at the end of each LFP+LMP 
convergence until the load flow failed Convergence was obtained for seven load mere 
ments When the system power level was incremented further, the load flow failed to 
converge The bus voltages and the corresponding real and reactive powers at the end 
of 7th load increment are given in Table 4 11 

The variation of the voltage and the P and Q loads for different states of the system 
are stored in files specified by the user The voltages, Pioad and Qioad of buses 12 and 
21 at different system loadings are shown in Table 4 12 




Table 4 10 Results at the end of LFP+LMP convergence 


Bus No 

Volt(p u) 

PJoad 

QJoad 

8 

1 0346 

6 0183 

2 0755 

9 

1 0488 

11 6264 

6 0339 

11 

1 0233 

7 7994 

2 8406 

12 

1 0068 

22 9119 

12 6432 

14 

1 0329 

6 4122 

2 4237 

15 

1 0302 

8 4535 

3 3825 

17 

1 0298 

9 1822 

5 5035 

21 

1 0255 

17 8055 

9 5469 

30 

1 0378 

10 8880 

3 7957 


Table 4 11 Results of the state prior to the failure of LFP 


Bus No 

Volt(p u) 

PJoad 

QJoad 

8 

9029 

12 5885 

4 4050 

9 

9240 

24 6422 

11 4605 

11 

9543 

17 0337 

6 5130 

12 

9513 

51 1413 

29 1964 

— 

14 

8881 

12 9331 

4 7384 

15 

8832 

17 0513 

6 6621 

17 

8880 

19 5865 

12 5561 

21 

8820 

37 9393 

16 2073 

30 

9047 

21 5883 

8 3055 
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T\ble 4 12 Voltage Pioad and Qio d at buses 12 and 21 different system loadings 


BUS 12 

A 

B 

C 

D 

E 

F 

1 

139 697 

44 445 

1 007 

22 912 

12 643 

2 

154 700 

49 615 

1 004 

25 493 

14 117 

3 

171 572 

55 434 

1 000 

28 432 

15 811 

4 

190 176 

61 610 

994 

31 768 

17 740 

5 

211 352 

68 712 

990 

35 657 

20 042 

6 

234 702 

76 180 

981 

40 089 

22 655 

7 

260 660 

84 083 

968 

45 190 

25 663 

8 

289 735 

92 636 

951 



51 141 


29 196 

BUS 21 

A 

B 

C 

D 

E 

F 

1 

139 697 

44 445 

1 026 

17 806 

9 547 

2 

154 700 

49 615 

1 019 

19 765 

10 454 

3 

171 572 

55 434 

1 011 

21985 

11 444 

4 

190 176 

61 610 

998 

24 444 

12 432 

5 

211 352 

68 712 

983 

27 262 

13 503 

6 

234 702 

76 180 

959 

30 394 

14 507 

7 

260 660 

84 083 

926 

33 914 

15 416 

8 

289 735 

92 636 

882 

37 939 

16 207 

(Tl 

iis data is stored by LOADMOD in the. file name 


specified by the user in file ’out_plot’) 

A 

System state number(first one is 

the base case) 

B 

Total P Joad of the system in MW 

C 

Total Q Joad of the system in MW 

D 

Voltage at the bus m p 

u 


E 

P Joad at the bus in MW 


F 

Q Joad at the bus in MW 
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The voltage of buses 4, 8, 11 and 12 with increased loading of the system are shown 
in Figiues 4 3 4 4 4 5 and 4 6 respectively Five conditions are taken to observe the 
system behaviour 

o Load modelling applied at 9 buses 

• Condition 1 + load mix change applied at bus 11 In this case the load mix 
change is applied at bus 11 in Models 6 and 28 Model 6 is commercial central 
an conditional , which consumes more reactive power whereas Model-28 is water 
heatei, which is a unity power factor load In each load increment the load of 
Model 6 IS increased by 0 1 MW and the same is deducted from that of Model 28 
Hence the voltage at this bus decrease as the total reactive power consumption 
IS increasing This effect can be seen in Figure 4 5 for the bus 11 In this case 
the decrease in voltage is small but when loads of all models (load mix) change, 
keeping the real power constant, the reactive power of all the models put together 
{Qioad at the bus) changes causing change in the voltage at the bus 

• Condition 1,2 + reactive power loss at the feeders taken into account and modelled 

• Condition 1,2,3 + synchronous condenser at bus 4 removed after 3 load increments 

• Condition 1,2, 3,4 + load at bus 12 is removed after 3 load increments 

Curves 1 to 5 in each of the Figures 4 3,44 45 and 4 6 represent, respectively, the 
voltages of the buses 4, 8, 11 and 12 for the above mentioned conditions The load was 
incremented in small steps as the system reached the stability hmit m order to capture 
the voltages properly 



Voingc at bus 4 



Total system real power load in MW 
ICEE 30 BUS SYSTEM 


Figure 4 3 Voltage at bus 4 for different conditions of loading 



Total system real power load m MW 


IEEE 30 BUS SYSTEM 


Figure 4 4 Voltage at bus 8 for different conditions of loading 
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Voltage at bus 1 1 



IEEE 30 BUS SYSTEM 


Figiue 4 5 Voltages at bus 11 for different conditions of loading 


Voltage nl bus 12 



Figure 4 6 Voltages at bus 12 for different conditions of loading 
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IEEE 30 BUS SYSTEM 

Figuie 4 7 Voltages at bus 19 foi different conditions of loading 

With the load modelling applied at nine buses it was observed that the voltages at 
buses 18,19 and 20 decreased rapidly when the loading was increased To capture the 
volt igos it tlicso throe buses accuiatcly they weie also modelled along with the nine 
buses The system loading was increased slowly to capture the static stability limit and 
it was obscived that the voltage at bus 19 was least among all Figuic 4 7 shows the 
voltage at the bus 19 for different loadings of the system Curve 1 shows the voltage at 
bus 19 without loads m the system modelled, and curve 2 is that with loads modelled 
at all 12 buses This shows that the actual system with load mix data as provided 
behaves as shown m curve 2 The voltage at bus 19 after considering the reactive losses 
and modelling them is as shown m curve 3 It is seen that the stability limit decreased 

Foi ihc same loidmg at the blis, the loactive powei calculated by the models foi 
diffeient load mix data is different because of which the modelled real and leactive 
losses also change This resulting reactive power causes the bus voltage to change 
when taken as input for the subsequent iteration of load flow Hence the load mix 
data IS the mam requirement for accuracy of final voltage levels With good load 
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St-xfcistics available when all load buses m the system are modelled, the results of load 
flow approach leahstic values and the ststic stability limit can be captured properly 



Chapter 5 


Conclusions 


• In the piesent work the characteristics of loads in the power system were consid 
eied for load flow and static stability studies The fixed compensation at the load 
bus was modelled for the changes m bus voltages Real and reactive losses in the 
feeders were modelled for the change in the loads for which the initial reactive 
power losses m the the feeders at 1 0 p u voltage at the terminal of the loads form 
the input 

• The results of load flow with load modelling depend on the load mix data For 
proper perfoimance, the data provided for the base case should be such that the 
real and reactive power outputs of all the models put together at 1 0 p u voltage 
plus the real and reactive power losses should be close to the specified P and Q at 
the bus If the load mix data is not proper then there is a imsmatch between the 
calculated and specified values of reactive power under base load condition This 
results in improper voltage levels at the buses Hence accurate load mix data is 
the mam requirement 

Strengths 

• Flexibility is provided for the user to study load flow with load modelling and 
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its application to static stability with different options for load modelling load 
mcicment iteration limits tolerance values for convergence of voltages generation 
of output files for each specified bus containing the information regarding the bus 
voltage real and reactive powers, etc 

• The load modelling part of the program is developed in terms of different modules 
and can be linked with any load flow program easily 

Limitations 

• The buses where only Z load is present (only compensation) and are required to 
be modelled for fixed compensation need 

— to be represented in P Q load data with zero P load and Q load in MW 
instead in Z load data 

- the specification of zero P load for all the models in the load mix data 

Hence m a big system having many Z loads, the load composition data of all the 
models at all these buses is required to be fed 

• The program is not supplied with any default load mix data 

Scope for further work 

• In the present work all load models are considered to be present at the bus and 
hence the votage input given to the models is the bus voltage However, voltage 
drops in the feeders may be modelled, and the resulting drops can be taken into 
into account for calculating the voltages at the models 

• The dynamically behaving equipments may he modelled for their extended use m 
transient studies 

• In countries like INDIA wliere the loads are high and voltage profile crosses the 
limits and coUapsts are frequent, load modelling with accnratc load mix applied 
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at all loid buses would help in accurately assessing the voltages and the margin 
of the stability limit this would require a survey to get the load mix and the load 
characteristics 
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Appendix A 

Computer code for Load Modelling 
Part of LOADMOD 


C itoH * 4oti >(> Xoio** % 4(4' IK 4 If 4' 4* ♦ 4i 4o»of >to|o^ If )|i iK % 41 left 
c LOAD MODELLING PART 

C OF 

C LOADMOD 

c MODULE 1 

Q iflf ifDflf If Iflf If If If If If If Ififlf If 4llf If 

c VARIABLES USED IN LHP 

real bpvr(lOOO) vOlClOOO) 10 v i dt vO qqClQOO) p q 
real pi ql p2 q2 v2 

real laot loadCEO) fact 1 xxl(lOOO) xx2(1000) xx3 

real ppload ppgon pploeaoa qqgen qqload qql asea absl abB2 abs3 

real Ira chong (1000 600) Im chan ab84 

integer bn(lOOO) tnb mmn ij 12 nail load itr nbli ii2 

integer 16 bnlic(600) load r q load itm(EO) tmpl tmp2 

integer load ralx y load mix itr Im buaea Im models 

integer Ira buano(lOOO) Ira mod no(BOO) Im marfcl Im mark2 

integer r pause 

oharactor*20 Ira char 

integer req out lile req all ol out unit 

integer out nobus out itr out buano(lOOO) 

oharactor*6 o 1 

oharacter*2 0 11 

integer req gen out no gen out gen out (1000) 
integer gen out Itr(lOOO) iS j8 
real gen out laot (1000) llagl23 
integer req Id out no„ld out Id out (1000) 
integer Id out itr (1000) Idobua(lOOO) marks 
real Id out lact(lOOO) Idolao(lOOO) Id outage pll 
real 12 111 112 rec buses rat bus (BOO) rec loss(BOO) 
integer req rec loss ii4 llagl lip Imp llplmp itr(lO) req vtol 
integer limit var lim Hag 
real v limit v lim temp(lO) v liraitl 
complex sKBOO) s2(600) 
c FILES USED IM - LHP - 

opan(unit BOBEO lile itr out! ) 

open(unit'«8989 lile* itr out2 ) 
open(xmit 6789 lile» out plot’) 
open(unit 6788 lile- out lile names ) 
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p n(unit 7000 file gen utag ) 
op nCunit 7001 fil load outag ) 
op n(unit 5010 file itr limit ) 


c FILE Ifdata CONTAINS THE INPUT DATA FOR lFP 
0 FILE nr out IS THE OUT PUT OF LFP WHERE THE BUS VOLTAGES 
c P & Q ARE PRINTED AT THE END OF EACH LOAD FLOW RUN 

c FILE data out IS THE OUT PUT OF LFP WHICH GIVES THE INPUT 
AND OUT PUT DETAILS 


c MODULE 2 

i(‘*+*>i<***+*****i|i**»** 


c 

0 qq IS A REAL ARRAY IN WHICH THE COMPENSATION AT THE LOAD BUS IS STORED 
c WHICH IS USED FOR ITS MODELLING IN LHP AND THEM SUBTRACTH) FROM THE 0 
qq(i) x3 


c MODULE 3 

If If ififif If ifif If iffifif if«if ifififf If 

c READING DATA FOR LIMITS IN ITERATIONS OF LFP LFP+LHP CONVERGENCE AND 
VOLTAGE TOLEREHCE FOR CONVERGENCE OF LFP+LMP FROM FILE tr limit 
readCBOlO ♦) 
roadCSOlO ♦)lfp 
road(6010 *) 
road(6010 *)lmp 
read(B010 ♦) 
r«ad(S010 ♦)v limit 
V limitl V limit 
roadCEOlO ♦) 
read(5D10 »)req vtol 
readCSOlO ♦) 
readCSOlO >*)limlt var 
roadCBOlO ♦) 
do 6011 1*1 limit var 
rendCBQlO *)lfplmp itr(i) v lira temp(i) 

SOU eontinua 

lim flag 1 

c READING DATA FOR GENERATOR OUTAGES DURING LOAD INCREMENTS AT THE 
c end of LFP+LHP CONVERGED RUN FROM FILE gen outages 
read(7000 ♦) 
roadCTQQO ♦)req g«n out 

ifCraq gen out aq l)tlian 
raad(7000 *) 
raad(7000 ♦)no gen out 
raad(r000 ♦) 
do 7200 i»i no gan out 

rBad<7000 ♦)out itr gen out(i) gen out itr(l) gen out f ct(i) 

7200 continue 

else 

clos#C7000) 

endif 

0 READING DATA FOR LOAD OUTAGES DURING LOAD INCREMENTS AT THE 

c END OF LFP+LMP CONVERGED RUN FROM FILE 1 ad outages 

r«ad(7001 +) 
raad(7001 ♦)req Id out 

if(req Id out eq i)then 
raad(7001 
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7100 


0 


c 


read(7001 *)iio Id out 
r ad(7001 *) 
do 7100 1 1 n Id 0 t 
readC700t *)out itr Id utCi) Id 
cont inu 
la 


ut itr(i) Id out lactCi) 


no Id out 0 


oloa (7001) 
ndif 

VARIABLE marks IS THE FLAG FOR NOTING THE LOAD 
DURING THE RUN TIME ITS VALUE CHANGES TO 1 VHEN 
OUTAGE OR ELSE IT REMAINS 0 


OUTAGE 
THERE IS AN 


marks 0 


c 


c 


c 


READING DATA FOR GENERATION OF OUT PUT FILES FOR V P 0 AT RIKPO 
FOR TOTAL REAL t REACTIVE POWERS OF THE SYSTEM THROUGHOUT THE LOAD 
INCREMENT ITERATIONS ( LFP+LHP CONVERGED ) FROM FILE out plot 


THIS USES ANOTHER FILE 
FILE NAMES ( o BuaNo 


ut fil name TO CREATE b STORE THE DEFAULT 
) WHEN THE ABOVE OPTION IS CHOOSEN 


read(6789 ♦) 
road(67a9 *) 
raad(6789 ♦)r8q ut Jil 


lf(r q out file «q Dthan 
r ad(6789 ♦) 
r<iad(6789 ♦)roq all of 

ifCraq all of oq l)than 
0 f 1 0 

do 655 1 1 nbua 
if(i It 10)than 
urita(6788 866)o fl i 

656 format(a2 il) 
endif 

if((i ga 10) and (i It 100))theii 
writ8(6788 6B8)o fl i 

658 format (a2 i2) 
endif 

ifCCi ga 100) and (i It 1000))th8n 
writa(6788 669)o fl i 

659 fortnat(a2 13) 
andif 

655 oontinue 

01088(6788) 

open(unlt 6788 fila out fila namaa ) 
out unit 6000 
do 657 1 1 nbua 

out unit out unlt+1 
rBad(6788 ♦)o f 
opan(unit out unit fila o f) 

657 continue 


ce>WTRAL Lr 

Msh A ^ ^ 


alae 

raad(6789 *) 
raad(6789 ♦)out nobua 
r8ad(6789 ♦) 
do 670 i 1 out nobua 
out unit 6000 

rand(6789 *)out itr out busno(l) o f 
out unit out unit+out buanod) 
opan(unit»out unit file o f) 

670 continue 


endif 


I 
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ndif 


c 


c 


WRITING THE 
writ (8989 
writ (8989 
wr te(8989 
writ (8989 
writ (8989 
writ (8989 
writ (8989 


OUT PUT DETAILS IN THE FILE 
♦)( ITN BUS NO VOLTAGE 

♦)( (KV) 

♦) 

) 

*)( +****♦**♦♦* loadhqd ITR 
♦) 

*) 


itr 


ut2 

P 

(MW) 


1 *** 


Q > 

(KVAR) ) 


♦ ) 


READING THE OPTION FDR THE CHANGE IN THE 
BUS FROM FILE 1 ad mix change 

p n(uait 9000 file lo d mix haage ) 
readOOOO ♦) 
r ad(9000 ♦)load mi y 


LOAD MIX OF THE LOAD AT THE 


if(l ad mix y e<j l)th n 
road(9000 ♦) 

road(9000 ♦)! ad mix itr 
elao 

cloae(9000) 

ondlf 


READING DATA FOR REACTIVE POWER LOSS IN THE FEEDERS AT THE LOAD 
0 BUSES FOR ITS MODELLING FROM FILE react 1 a 
opan(unit 321 file react 1 a ) 
raad(321 *) 

raid(321 »)req rec loaa 

if(req rec loaa eq Dthen 
r ad(321 ) 

read(321 ♦)rec buaea 
r ad(321 ♦) 
do 322 1 1 rec buaea 
read(321 *)rec bua(i) rec 1 s(i) 

322 contlnu 

elaa 

oloBe(321) 

endif 


c THE LOAD AT A BUS IS DIVIDED AMONG THE DEFINED LOAD MODELS 

c + r OSSES THE E DETAILS OF ALL BUSES IS IN FILE load compoaitiou 

opea(unit 31 file load compoaitiou ) 
uiirm 31 

c THE FILE load modelling IS THE MAIN OPTION FOR LOAD MODELLING 
0 WHICH ACTIVATES LHP - THE DATA FOR CHOICE OF BUSES AND 

c THEIR REAL POWER LOADS ARE READ FROM THIS FILE AND LHP IS 

c ACTIVATED TO GET THE REAL AND REACTIVE POWER AT THESE BUSES 

c WHICH ARE STORED IN ARRAY al THESE ARE USED FOR MODELLING 

0 REAL AND REACTIVE POWER LOSS IN THE FEEDERS AT THE LOAD BUSES 
open(unit 5 file load modelling ) 
road(5 ♦) 
r ad(5 *)ttab 
read(S *) 
raad(B ♦)fO f dt 
raad(B *) 

if (tub eq 0)then 
rewiud(E) 
elae 

do 7991 1 1 tub 
read(B ♦)bn(i) bpwr(i) vOl(i) 

V"1 0 

1 SO 0 
dt 0 0 
vO V 
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fO 1 
ij bn(i) 
p 0 0 
q 0 0 

f ct 1 1 0 
Im mark2 0 
Im mod 1 1 

Im mod no 1 
Im change 0 0 
m rk3 0 
Id outage 1 0 
sl(ij) (0000) 
flagl 1 

call modelCv f dt vO fO mmn i ij p q fa t 1 Im mark2 
* Im m d la Im mod no Im hang marks Id utair 1 fl bI 

Klj) cmplxCpll q) ® ^ ’ 

7991 c ntinuo 

cl a (5) 
ondif 

cloao(nnnn) 

fact 110 
flagl 0 

c WHnN LFP+LHP CONVERGES AND LOAD IS INCREMENTED FOR NEXT STATE 
c OF THE SYSTEM IT RUNS FROM THE STATEMENT LABELLED 8980 
c VARIABLE l6 KEEPS TRACK OF THE ITERATION NUMBER 

c WHEN ONE LFP+LMP RUN IS OVER IT IS REPEATED FOR NEXT RUN FROM 
c STATEMENT LABELLED 899 VARIABLE i2 KEEPS TRACK OF THE 
0 ITERATION NUMBER 
16 0 

8980 12 0 

16 16+1 

899 12 12+1 

write (8989 ♦)( (LFP + LHP) ITR i2 ) 


c MODULE 4 

c THE FILE WITH UNIT NUMBER 2 ( nr out ) IS CLOSED AT THE END 

c OF EACH LOAD FLOW RUN 
cloae(2) 

c THE DETAILS OF THE BUSES TO BE MODELLED ARE READ FROM 1 d modelling 
c AND LOAD COMPOSITION FROM load oompoaltion 
open (unit B file load mod lling ) 
read(B *) 
read(B ♦)tnb 
raad(B ») 
read(S ♦)fO f dt 
r9ad(B *) 

0 +++++++++++++++++ 
if(tnb *q 0)then 
rewind(S) 
elae 

do 799 1 1 tnb 

r9ad(B ♦)bn(i) bpwr(i) vOl(i) 

799 continue 

cloae(B) 

open(unit»31 fil load_c mpoaition ) 
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fO f 
ij bn(i) 

p 0 0 

q 0 0 

fact 110 
Im ma k2 0 
Im mod 1 1 

Im ro d n 1 
Im Chang 0 0 
marks 0 
Id utag 1 0 
sl(ij) (0000) 
fl gl 1 

call modalCv f dt vO fO nmm i ij p q f ct 1 Im marks 
* Im m d Is Im mod n Im chang marks Id out g 1 fl gl pll) 

sl(ij) mplxCpll q) 

7991 continu 

closoCS) 
eridif 

close (nnun) 

fact 110 
flagl 0 

c WHEN trP+LHP CONVERGES AND LOAD IS INCREMENTED FOR NEXT STATE 
c OF THE SYSTEM IT RUNS FROM THE STATEMENT LABELLED 8980 
c VARIABLE l6 KEEPS TRACK OF THE ITERATION NUMBER 

c WHEN ONE LFP+LMP RON IS OVER IT IS REPEATED FOR NEXT RUN FROM 

e STATEMENT LABELLED 899 VARIABLE i2 KEEPS TRACK OF THE 
0 ITERATION NUMBER 

l6 0 

8980 i2”0 

16 i6tl 

899 i2 i2+i 

writa(8989 ♦)( (LFP + LHP) ITR i2 ) 


♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦it!* 

MODULE 4 

: THE FILF WITH UNIT NUMBER 2 ( nr out ) IS CLOSED AT THE END 

: OF EACH LOAD FLOW RUN 

cloae(2) 

c THE DETAILS OF THE BUSES TO BE MODELLED ARE READ FROM 1 d m d llmg 
c AND LOAD COMPOSITION FROM load composition 
openCunit B file load mod lling ) 
readCS ♦) 
readCS *)tnb 
readCS *) 

X ad(S ♦)fO f dt 
read(B *) 

C •f++4"t'+-t'+++++++4"f+ 

if (tab eq 0)then 
rewiad(8) 

else 

do 799 i 1 tub 

readCB *)bn(i) bpwr(i) vOl(i) 

799 continue 

cloae(S) 

openCunit 31 file load_composition ) 
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nimn 31 

pi 0 0 
p2 0 0 

ql 0 0 
q2 0 0 


c 


488 

459 


READING THE DATA FOR THE CHANGE IN THE LOAD HIX 
BUS FROM FILE 1 ad mix change 
if(12 eq l)th n 
if (load mix y eq l)then 
r ad(9000 *) 
r ad(9000 *) 

r ad(9000 ♦)lm char Im buses 
r ad(9000 ♦) 


OF THE LOAD AT THE 


do 459 1 1 Im buses 
read(9000 *)lm char Im buano(i) 
r ad(9000 ♦)lm char Im m d la 
do 458 j 1 Im mod Is 

r0nd(9OOO ♦)lra char Im mod no(j) Im har Im chan 
tmpl Im busnoCi) 
tmp2 Im mod 110(3) 

Im Chang (tmpl tmp2) Im chang (tmpl tmp2)+lm han 
continue 
continue 
ondif 
endif 


Im markl 1 


0 EACH BUS TO BE MODELLED IS TAKEN ONE BY ONE 

e BY TAKING INTO ACCOUNT ALL THE OPTIONS 

do 798 i 1 tnb 

ll((load mix y oq 1) and (Im busnoClm markl) eq bn(i)))then 
c VARIABLE Im marlE2 IS THE FLAG FOR CHANGES IN LOAD MIX ITS 

0 VALUE IS MADE 1 IF LOAD MIX IS PRESENT OR ELSE REMAINS 0 

Im m rk2 1 

if (la markl It Im bu3es)thon 
Im markl Im markl+1 
andif 
else 

Im mark2 0 
andif 

vO vOi(i) 
ij bn(i) 

c TRANSFERING THE VOLAGE AT BUS ij OBTAIHD FORM LOAD FLOW PART IN 

0 VARIABLE X TO VARIABLE v TO BE USED IN LOAD MODELLING PART 

V x(nbU8+ij) 

marks 0 

if(roq Id out eq l)then 
do 7030 i8 1 no Id out 
if((i6 1) ga Id out itr(i8))than 
j8 Id out(i8) 
il(j8 eq iJ)thon 
Id outage Id out fact(i8) 
maxk3 1 

pload(j8) pload(j8)*ld outage 
qload(j8) qload(j8)*ld outage 
go to 7031 
else 
marks 0 
Id outage 1 0 
andif 
else 
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marks 0 

Id outag 1 0 

ndif 

7030 ntinu 
also 

marks 0 
Id out g 10 
BUdlf 

7031 flagl23“0 0 


2 V 

V v*v0 


FOR EACH BUS SUBROUTINE mod 1 IS CALLED WHICH CALCULATES 
FINAL P t q OF THE LOADS AT THE BUS 
call modalCv f dt vO fO tinnn i ij p q fact 1 Im mark2 
* Im models Im m d no Im hanga marks Id outage si flagl pll) 


if (marks q Dthan 
qq(ij) 0 0 
endif 

PUTTING BACK CHANGED REAL ft REACTIVE POWERS AT THE BOS OBTAINED AFTER 
! APPLYING THE LOAD MODELLING FOR CHANGED VOLTAGE INTO THE VARIABLES 
: pi ad ft ql ad WHICH ARE USED BY LOADFLOW PART 

pload(ij) p 

: MODELLING THE COMPENSATION 

qqClj) qq(ij)tv2>fv2 

qload(lj) q qq(ij) 

MODELLING THE REACTIVE POWER LOSS IN THE FEEDER AT THE BUS 
if(bpvr(i) gt 0 OOOODthan 
ifCreq rao loss eq l)than 
do 324 ii2 1 rec buses 

lf((rac bus(ii2) eq ij) and (marks q 0))tlien 
a2(ij) iiipli(pll q) 

112-cabB(B2(i3))/caba(sl(ij)) 

111 r c lo 3(ii2)/(v2*v2) 

12 Illt.ll2tll2 

qlond(lj)»qlond(ij)+12 

endif 

324 continue 

endif 
endif 

marks 0 

c WRITING THE MODELLED BUS VOLTAGE P ft q IN OUTPUT FILE itr out2 
write(8989 891)12 bn(i) v pload(ij) qload(i3) 

891 formatd* iS 3x i3 7x f8 4 3x f9 4 3x f9 4) 

0 VARIABLES p«l t qinj ARE INJECTIONS OF REAL ft REACTIVE POWERS 

0 AT THE BUS USED IN THE LOAD FLOW PART 

pinj(ij) pload(ij)/baBe 
qlnj(ij) qload(ij)/baae 

798 continue 

cloBe(muui) 

oloBoCnin) 

c STORING PERSENT VOLTAGES IN ARRAY xxl FOR VOLTAGE TOLERENCE 
c COMPARISIQN M THE NEXT ITERATION 
if(i2 eq Dthen 
do SOSO i“l abUB 

XXlCi)“x(ttbnB+l) 
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5050 c ntiau 
endif 

STORING THE MAXIMUM VALUE DIFFERENCE IN THE VOLTAGE IN VARIABLE xx3 
xx3 0 0 

ifCa.2 go 2)thon 
d 50S1 1 1 nbu 

xx2(i) ab ( xxlCi) x(nbu +i) ) 
if(xx2(i) gt xx3)th n 
xx3 xx2(i) 
sndif 

5051 c ntinu 

do 5052 1 1 nbua 
xxl(i) x(nbu3+i) 

5052 continuo 
ondif 

c SETTING THE VOLTAGE TOLFRENCE LIMIT IF IT CHANGES 
lf(raq vtol aq l)thon 
if(i2 oq Ifpltap itrClim flag))th u 
V limit V lim t mp(lim flag) 
lim fl g lim flag+1 
andif 
ondif 

if(i2 gt Imp) then 

writoC* ♦)( COULD NOT CONVERGE IN Imp ITERATIONS OF 
* LOADHODEL PROGEAM ) 

atop 
ondif 

write (♦ *)( END OF LOAD MODEL t LOAD FLOW ITR i2) 
writo<* ♦) 

write (* *)( VOLTAGE DIFFERENCE IN ITERATIONS IN P U xx3) 
writo(f ♦)( VOLTAGE LIMIT TOLERENCE IN P U v limit) 

write (♦ *) 

if((xx3 go V limit) or <i2 oq l))thon 
go to 899 
elao 

write (♦ ♦)( ^ 

write(*»)( LFP + LMP CONVERGED) 

wxitaC* ♦)< ^ 

ondif 

endif 

C ++++++ 

lim flag 1 
V limit V limit! 

write (8989 *) , ^ 

write (8989 *)( ♦*+*+++*♦♦ LODMOD ITR i6+l *********** ) 

writo(8989 ♦) 


0 

c 


0 




WRITING THE VOLTAGES REAL AND 
IN THE OUT PUT FILE itr outl 
THE CONVERGENCE OF LFP+LMP 


reactive power injections 
THIS IS WRITTEN AFTER 


AT ALL BUSES 


ppgen 0 0 
ppload 0 0 
qqgetl 0 0 
qqload 0 0 
pploeaos 0 0 
qqloBS 8 0 0 


do 60561 i 1 nbus 
if(iflag(i) eq Dthon 
ppload ppload+pinj(i) 
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qqload qqlo*id+qinj (i) 

Is 

PPg n PPg v+pinj(i) 
qqgon qqgen+qinj (i) 
endlf 

pploaaoa ppg n+ppload 
qqloaa qqg n+qqlo d 
B0551 continua 

WRITING INTO GUT PUT FILE itr outl 
writeCSOSSO ♦)( ITR No i6 

writ (S0B60 ♦)( Baa KVA ba e Pg 
wr te(S0550 *)( Pload ppload Los 
writ (50S50 *)( Bus No V Pin] 
write (60S60 *)( (p u) (p n) 

d SOBS 1 1 nbus 
writoCBOSBO B0BB2)i i(nbus+i) pinj (i) qin](i) 
BOSB continue 

B0BB2 Iormat(i3 2x 3(f9 4 2x)) 


) 

PPgen) 

PPl ) 
Q1 ad ) 

(p u) ) 


c WRITING THE VOLTAGES AND P A q INJECTIONS AT THE BUSES IN SEPERATE 

c FILES AS PER THE REqUIREHENTS SPECIFIED IN DATA FILE ut pi t 

lf(raq out lila eq l)than 
llCraq all of oq l)th n 
out unit 6000 
do 681 1 1 nbus 
out unit out xmit+1 
abal aba (ppload) *baaa 
abaq aba(pinj (i))tba3a 
abs3 aba(qinj Ci))*base 
abs4 aba (qqload) ♦base 

WriteCout unit 682)i6 abal aba4 x(nbu +i) abs2 absS 

681 continue 

682 format(i3 B(3x f9 3)) 
else 

do 683 i 1 ut n bus 
out unit 6000 

out unit out unit+out bu no(i) 

] out buano(i) 
abal aba (ppload) ♦base 
abs2 abB(pinj (j))^baaB 
abs3 ab8(qinj (j))^base 
aba4 aba(qqload)*base 

write (out unit 684)i6 absl aba4 x(nbu +]) abs2 abs3 

683 continue 

684 formatdS B(3x f9 3)) 
ndif 

endif 

c GENERATOR OUTAGES 

if(req gen out eq Dthan 
do 7010 i 1 no gen out 
if(i6 ge gen out itr(i))then 
jagen out(i) 

pgett(j) pgan(j)*gan out fact(i) 
qmax(j) qmax(j)+gen out_fact(i) 
qmin(j) qmin(j)*gan out fact(i) 

C WHEN GENERATORS ARE OUT THEY ARE TREATED AS LOAD BUSES 
c FROM THE NEXT ITERATION BY MAKING IFLAG(BUS) ZERO 

0 IFLAG IS A VARIABLE IN LFP WHICH IS 

e 1 FOR LOAD BUSES - 

c 2 FOR GENERATOR BUSES 

c 3 FOR SLACK BUS 

iflag(j) 1 
c 


endif 
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7010 continua 

and f 

c LOAD OUTAGES 

ifCroq Id out aq l)then 
do 7020 1 1 no Id out 
if (16 go Id out itr(i))then 
j Id out(i) 

ploadCj) pi d(])H<ld out fact(i) 
ql ad(j) qload(j)*ld out fact(i) 

Id bus(i) 0 

Id fa (l) Id out fa t(i) 
marks 1 
ndif 

7020 nt nu 

endif 

c READING THE DATA FOR LOAD INCREMENT FROM FILE 1 d in ram nt 
if<i6 q Dthen 

openCunit 606 file load incr mont ) 
raad(50B ♦) 
roadCSOB *)load raq 
if (load raq aq 0)thon 
atop 
endif 

r adCBOS *) 
raadCBOB *)load itr 
raadCSOB *) 
raad(B0S *)r pau 
raad(50B ♦) 
do 609 1 1 load itr 
raadCEOS *)load itm(i) fa t 1 ad(i) 

609 continua 

read(506 ♦) 
raad(B0S ♦)nnll 
raadCBOB *) 
r ad(S0B ♦)nbli 
raadCSOB ♦) 
do BOS 1 1 nbli 
raad(B0B *)bnlic(i) 

608 continua 

endif 

0 multiplying the load FACTOR WITH THE NEW AHHOUNT OF INCREMENT 
fact 1 fact l*faot loadCiB) 


c 


0 

c 


326 

606 


EASING THE GENERATION AND LOAD 
nail q Dthan 
I E06 i 1 nbua 

igen(i) pgan(i)*fact load(16) 
)load(i)»>pload(i)'t'fact loadCiS) 
lload(i) qloadCi)H‘faot loadCiS) 
qmaxCi) qmax(i)*fact lond(i6) 
qmin(i) qmin(i)*f ct„load(iG) 
if (raq rec loan aq l)than 
do 325 ii4 1 rac buaea 
if(rac btt8(ii4) aq i)tben 
rec loa»(li4) rec lo»fl(ii4)*fact loaddS) 

endif 

continue 

endif 

continue 


B07 i»l nbli 


j-bnUc(i) 

pgen(j)”pgon(j)’t'fact load(i6; 
pload(j)'»pload(j)*fact load(i6) 

qload(j) qload(j)*faot load(i6) 
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qm x(]) qm x(3)*fact Ioad(i6) 
c qmin^j) qmin(j)ffact lo d(i6) 

S07 continue 

ndif 


ii(i6 It 1 ad itr)than 
«rito(* *) 

Hrita(* )< 
writo(» ♦)( 

Mrito(* *)( 
vrite(* ♦)( 
writaC* *)( 
if(r pans 
pause 
ondif 

go to 8980 
endif 


* , 

* ) 

END OF LOAD IHC ITR 
* ♦ ) 
l)thBn 


16) 


writaC ♦)( 
writaC* +)( 
«rlte(* ♦)( 
writoC* *)( 
writaC* *)( 


♦*************»******^*^^^^^^^ ^ 

* * ) 

END OF LOAD INC ITR i6) 

* ♦ ) 


ifCr pause aq Dthon 
pause 
andlf 
closaCSOS) 


*! 




c SUBROUTINE FOR READING THE MODEL LOADS 

c IKPLEHENTINC INCREMENTS OUTAGES AND MODELLING THE REAL LOSSES 

subroutine model Cv f dt vO fO nunn iji ijij p q fact 1 

* Im iQark2 Im models Im mod no Im chang marks Id utage 

♦ al flagl pll) 

implicit none 

integer no model iji ijij ImmS Imlm marks flagl 

integer i an naan Im marks Im models Im mod noCSOO) 

real realp(lOOO) reaotp(lOOO) Iw ehangeClOOO BOO) 

real V vO f fO dt add change 113 114 

real p q cat loss mlClOOO) fact 1 Id outage v21 pU 

character*20 area cl c2 

complex slCBOO) sSCEOO) 

c ~ 

0 opstt{unlt«30 file* program dat ) 

0 openCunit 32 file program data ) 

open(unit**44 file* progtaml out4 ) 
opan(unit 4S file programl outS ) 
open(ttnlt«22 file* comp code ) 
c _ _ 

0 no model gives the total nuaer of models used 

cat loss 0 0 

0 READING INFORMATION REGARDING NUMBER OF MODELS FROM FILE comp code 
read(22 *) 
read(22 ♦) 
read(22 *)no„modol 
reHlnd(22) 

0 READING THE NAME OF THE AREA FROM FILE load composition 
readCmmn *) 
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r adCiuinn 13)aroa 
13 format Ct38 a20) 

raadCnimn *) 

lmm3 1 

do 1011 i 1 ao modal 

READING THE MODEL LOAD INTO THE VARIABLE ml 
roadCnnim *)nn ml(i) 

C NOTING LOAD MIX CHANGES 

If (dm raark2 aq 1) and (Im mod no(lnim3) q i))th n 
Imlm Im mod no(lmm3) 

add Chong Im chang (ijij imlm) 
lf(lnim3 It Im models) then 
lmm3 lmm3+l 
ndif 
else 

add chonga 0 0 
ndif 

c IMPLl MI N riNC MQDEI LOAD INCREMENT IN LFP+LMP CONVERGED ITERATION 
c LOAD MIX CHANGE AND LOAD OUTAGES 
mlCx) ml(i)*fact 1 
mlCi) ml(i)+add change 
mlCi) ml(i)*ld outage 

1011 continue 

c READING THE LOSS AND ITS HODEFICATION FOR LOAD OUTAGE 
road(nnnn *)ol cat loss c2 
if (marks oq l)th n 
cat loss cat laaB*ld outage 
andif 

ifCiji aq Dthan 
ifCflagl «q 0)then 

wiita(45 *)( THE VARIATION IN P 4 q WITH VOLTAGE 4 FREQUENCY IN ) 
168 format (20x 20a) 

writa(4B *)( BUS TIME VOLT FREq P q AREA ) 

ondif 

andif 


c THE MODEL LOAD INFORMATION PRESENT VOLTAGE FREQUENCY TIME STEP 

0 ARF TRANSFERED TO SUBROUTIN modelehar WHICH IMPLEMENTS EACH MODa 

C CHARACTl RISTICS THE REAL AND REACTIVE POWERS OF ALL MODELS AFTER 
0 THE CHARACTERISTICS APPLIED ARE STORED IN REAL ARRAYS r alp 
0 AND raaotp 

call modeloharCv vO f fO realp reactp dt ml) 
ifCflagl eq 0)then 

c write (44 *)(. THE LOADS OF MODELS WITH CHARACERISTICS ) 


c 

write (44 169) area 
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format (2x AREA 

Sx 20a) 


c 

write(44 •)( Voltage 

V Frequ ney 

f) 

c 

write (44 



c 

writa(44 ♦)( MODEL NUMBER REAL POWER 

REACTIVE POWER ) 

c 

write (44 ♦)( fact 1 

fact 1) 


<J 

write (44 ♦)( 

( in MW ) 

( in MW ) ) 

c 

write (44 *) 




andif 


poo 

q 0 0 

c ADDING ALL THE MODELLED REAL AND REACTIVE POWER OUT PUTS INTO 
= VARIABLES p and 'q 
do 172 i 1 no model 



ii(fl gl oq 0)th n 

c writ (44 93)i r alp(i) r actp(i) 

endif 

p p+roalpCi) 
q q+r actp(i) 

172 continue 

pll p 

93 format(5x i3 tl5 fll 4 t31 ill 4) 

C MODELLING THE LOSSES AT THE FEEDERS 
ifCp gt 0 OOOODth. n 
ifCflagl oq 0)than 
v21 v/vO 

aSCijij) cmplxCp q) 

114 cabaCa3(iji3))/cab3( l(ijij)) 

113 cat 1 s /(v21* 21) 
c t loss 114*114*113 
1 a 

oat loaa 0 0 

ondiT 

andii 

ifCflagl q 0)then 
0 writ (44 ♦) 

c wrlto(44 99)cat loaa 

99 f rmat(2x LOS ES 3x flO 4 HH ) 

endif 

0 LOSSES BEING ADDED TO TOTAL LOAD OBTAINED AT THE BUS 
p p+eat loaa 

ifCflagl eq 0)then 
c write (44 ) 

c write (44 ♦)( TOTAL REAL POWER p MW REACTIVE POWER 

c * q MW ) 

0 write (44 *) 

0 writo(44 *) ( « ) 

c wxa.ta(44 *) 

endif 

ifCflagl eq 0)then 
wrlt«(4B 166)ljij v f p q area 
166 format (13 lx f9 3 lx f6 3 lx 19 4 lx 19 4 lx alO) 

endif 

return 

end 

*♦»*♦♦♦***♦*♦♦*♦**♦*•**♦♦*♦♦♦*»♦♦♦****♦*♦♦*♦♦*♦♦+♦**♦♦♦* ♦ 

c SUBROUTINE FOR 

c IMPLEMENTING THE MODEL CHARACTERISTICS 

+>♦<*# eeeeeeee+eeeee+eeeeeee*# ★♦'♦fete #*********♦**♦+♦♦♦*** ******* 


aubrovitine modeloharCv vO f 10 roalp raa tp dt ml) 
implicit none 
Integer a d 1 ii 

integer nroodel model no nom nan flag nlile 

real theta realp(lOOO) reactp(lOOO) 

real p q pO qO v vO 1 10 dt v2 ml(lOOO) 

real model load modn nr rm nmm dp dq 

real pol pv pf qv qf n« pofnm pvnm plnm qvmu qfnm 

character*! e p exp pol 

character*6 model code m c 

openCunit 41 file a d ) 
openCunit 42 file comp char ) 



p nCunit 49 filo comp harp ) 
openCunit 51 filo d d ) 
nfilo 51 

r ad(4l *) 
r d(41 *) 
roadCdl f)!!!!! dal 
road(41 ♦) 
r ad(41 *) 

roadC42 ♦) 
road(42 *) 
read(42 *) 
r ad(42 *) 
raadC42 *) 

nn 0 0 

EACH MODEL IS TAKEN AND THE OPTIONS FOR THE TYPE OF CHARACTERISTICS 
EXPONENTIAL OR POLAR IS READ FROM FILE a d 

do 171 i 1 nmodol 

read(41 ♦)fflodal no a d a p 
modal load 
flag 0 

« * «ii' * *♦»♦♦♦♦♦ 

if (a gt 0)than 
flag 1 
axp a 
pol p 

if (a p aq a3cp)thon 

READING THE PARAMETERS OF THE LOAD MODEL FOR EXPONENTIAL MODEL 
FROM FILE comp char 

raad(42 ♦)nom modal coda pof pv pf qv qf nm 
♦ pofnm pvjun pfnm qvnm qftrni 

lf(nom aq i)than 
raad(49 *) 
elaa 

vritaCf ♦) ( reading from the filea la wrong ) 
endif 

OBTAINING THE INITIAL Qo OF THE MODEL FROM IT S POWER FACTOR 

p 0 0 
<1 0 0 

pO modal load*(a»0 01) 
theta acosCpof) 
qO p0*tan(thata) 

IMPLEMENTING MODEL CHARACTERISTICS 

ifCCnm aq 0 0) or (nm aq 1 0))then 
For Nm 0 0 or Nm 1 0 

p pO*((v/vO)**pv)'t'((f/fO)**pf ) 
q qO*C(v/70)**qv)*(Cf/fO)**qf) 

alaaifCCnm gt 0 0) and (nm It 1 0))th n 
For Nm hetwaan 0 and 1 

p » nm*pO*((v/vO)**pv)#((X/iO)*apf) + 

♦ (1 nm)YpO*((v/vO)**pvnm)t((f/fO)’'"'‘P^““) 
q •• n«*qO+((v/vO)n<H‘qv)*((f/XO)**qi) + 

* (1 nm)aqO*((v/vO)**qvnm)*((f/XO)*»qfna) 

endif 

raalp(i)"p 
reactp(i) q 


alsaif(e p aq poDthan 
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C 

raadC49 *)n m and th oth r things 

c iJ(nom oq i)than 

c roadC42 ♦) 

c ol 0 

writ (♦ ♦)( reading from the file i wr ng ) 
ndif 

c THE POLYNOMIAL PART CAN BE ADDED HERE 

endif 

c 

endif 

c a =s 

ifCd gt 0)then 

c THE DYNAMIC PART CAN BE ADDED HERE 

if (flag oq O)thon 
read(4? ♦) 
road(49 *) 
ondif 

p 0 0 

q 0 0 tit 

0 READING DATA FOR DYNAMIC MODELLING OF INDUCTION MOTOR FROM FILE d d j 

readCBl ♦) 

readCSl *) ** 

readCBl *)m c m dn nr 

rand(6i ♦) l] 

v2«v/vO I 

do 1001 ii 1 nr 
road(61 ♦) 
r ad(Bl ♦) 

read(Sl *)rm nwm nnn 
read(61 *) 
xm 0 001tirm*nmm 
rm rra*Cd*0 01) 

call dynmodal(v2 f dt nfile dp dq) 

0 converting into MW 
p p+Crm*dp) 
q q+(rm*dq) 

1001 continue 

roalp(i) p 
reaotpCi) q 

endif 

c » »• 

flag 0 

171 continue 

c1obo(41) 
olosaC42) 
closo(49) 
oloso(Bl) 

return 

end 



Appendix B 


Dynamic modelling of induction 
motor 


In the dyn nine moddlmg of induction motor, the clasical model is taken and the swing 
equation is solved The equivalent circuit of induction motor is shown in Figure B 1 


R-s 

AW- 


Xs 


1 Xr 


Vs 


Xm 


Rr 

s 


Figure B 1 Equivalent circuit of induction motor 


After getting the Thevinin equivalent of the left of line 1-2 we have the circmt as shown 
in Figure B 2 This is true only if the supply voltage is strong, i e , V', is an ideal voltage 
source 
Here, 


jXmVs 

^ Ra+jiXu+Xm) 


(B 1) 
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Ri- Xe X 

AAA/' — '''757575''^ ^757575~\ 



I 


rigui6 B 2 EquivalGnt circuit of induction motor with. Thevenin equivnlent 


H jXg — 


J^Tn.{Ra + jXj) 

Ra+j{X +XJ 


Ir = 


K 

{Re + +j{Xe + Xj.) 


The thioc plnsc clectiical torque, 


o r 2 




where Pf — number of poles m the motor 
Substituting the value of /, m the above equation 


(B 2) 

(B3) 


(B4) 


T, 


= ft ) 

2 ^'■u>.-u,JiR, + ^r + {x. + x,y 


(B 5) 


ConvGitmg the angular fiequcncies u, and tOm into mechanical units, i e 


w”* =: uja/{Pf/2) and we get 


3i?r ^ 

i^T~UJ^\R.^^fP{Xe-VXrf 


(B 6) 


The swing equation for the motor is 


dt 


1 

2H 


{Te~Tm) 


(B7) 
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The mcchmic xl toiquc 

T,n = Tmo ( A«)2 + Bu^ + C) (B 8) 

wlieie A, B xnd C iie constants and are different for different motors 
Substituting the values of Te and m the equation B 7, we get, 

rf< ^ 1 ,, 3fi, Yl 

dt 2HK:--,^s(R, + :^r + lx, + X.y 

-r„„{yl«)^ + + C)] (B 9) « 

I 

Equation B 9 is solved using Runge Kutta fourth order method « 



Appendix C 


Computer code for the dynamic 
model of induction motor 


♦****»♦*♦♦♦♦*****♦**»*♦*♦»♦* ♦**)(l******+***»»**,|,******;^**,^*^*^^^^^ 

SUBROUTINE FOR 

0 DYNAMIC MODEL OF INDUCTION MOTOR 

subroutin dynm d«l(dva dl dt nfilo dp dq) 

implicit none 

complex aa bb co dd 

integer i niile 

real dvfl dt dt vo 

real dp dpfrac dsl dpf 

real rs xa xm rr xr re xe pi 

real dli da db dc 

real vba wa wm wmO tmO tm wmOi wm2 ffl 
real kC4) al a2 a3 a4i a42 a4 aS dq thetai 

tm 0 0 
wm 0 0 

c REAOIND DATA FOR INDUCTION MOTOR FROM FILE d d 

xeadCntile *)dp dpfrac daf dpf 
read(nfile *)ta xa xm rr xr 
read(nfilB ♦)dh da db do 

pi-22 0/7 0 

wba 2 0»pi*50 0 
wa 2 0+pi+df 

c - wmOl IS THE VALUE OF wm AT FULL LOAD AT i 0 p u 

0 VOLTAGE AND 60 0 Hz FREQ 

wmOl 2 O+piaSO Ofd 0 daf) 

c wmO IS THE INITIAL VALUE OF wm FOR THIS ITERATION 

0 - FROM NEXT ITERATION ONWARDS THIS INITIAL VALUE IS THE 

0 VALUE OF wm OBTAINED AT THE END OF PREVIOUS ITERATION 

0 - WHICH IS REQUIRED TO BE IMPLEMENTED 

wmO wmOl 

0 - CONVERTING INTO BASE VALDES 

wa wa/wba 
wmO wmO/wba 
wmOl wmOl/wba 
wm-wmO 
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dp dp*dpfric 
tmO dp/wmOl 

finding THEVENIN EqUIVALENT OF STATOR + MAGNETIC CIRCUIT 

a (0 0 1 0) 
bb=aa*xm*(r +(aat‘xa)) 
cc-rs aa*(xa+xm) 
dd-bb/cc 
r wroal(dd) 
xa airaag(dd) 

TUEVENIN EQUIVALENT OF VOLTAGE 
va ( a*xm*dva)/(rs+aa+(x3+xm)) 

- APPLYING RANGE KUTTA 4th ORDER METHOD 
do 100 i 1 4 
ffl 0 0 

al 1 0/(2 0*dh) 
a2 (3 Q*rr)/(wa win) 
a3 va*va 

a41 (rB+(rr'i'wa)/(wa wm)) 
a42 xa+xr 

al (a41*a41)+(a42*a42) 

wm2 1 0 (wm/wmOl) 
aE (da't'wm 2 *wm 2 ) + (db*wm2)+dc 

lit al’t'(((a2*a3)/a4) (tm0*aS)) 


k(i) dtffJll 


100 

c 


if ((i «q i) or (i «<1 2 ))thett 
wm tfm+(k(i)/2 0) 
alflaifd aq 3)thaii 
win wm+k(i) 

endif 

oohtinua 

RESTORING THE VALUES OF wm AND tm 
wm wmO+Cl 0/6 0)*( k(l) + 2 0*k(2) + 2 OfkO) + k(4) ) 
wm2 1 0 (wm/wmOl) 

tm tmO*C da*wm2*wm2 + db*wm2 + dc ) 

dpotmfwm 

daf (w8 wffl)/wa 

wmO wm 

thatal acoa(dpf) 
dq dp*tan(th8tal) 


taturn 

and 



Appendix D 

Data and result files 


File comp coda 


1 

Total number of modalla 

28 

agriup AGRIcultural Water Pump 

2 

arclur ACR FORaaea 

3 

centao CENTral Air C nditioning 

4 

cltdry CLoTbes DRYer 

5 

cltuer CLoTha WaSbeR 

6 

cmceao CoHmartial CEntral Air Conditi ning 

7 

ciurmac CoHmartial RooM Air Conditioning 

8 

colrtv CQLouR TalaViaion 

9 

oommhp CQMMartial Haat Pump 

10 

dabwar DISH WaSbaR 

11 

altrya ELacTRolYaiS 

12 

llrlit FLuoRaacent LIgbTing 

13 

fraz r FBEaZER 

14 

fuxfan FUEnac FAN 

15 

bpcaac Heat Pump CEntral Air Conditioning 

16 

bpcmao Haat Pump CoHmartial Air Conditioning 

17 

bpapbt Haat Pump SPace HaaTing 



18 

inolit INCandeacant LIgbTing 

19 

Iridlm LaRga InDuatriaL Motor 
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20 

ordfan 

ORDinary FAN 

21 

ordmtr 

ORDinary M T R 

22 

rdpmp 

ORDinary PuMP 

23 

pplaux Power PLant AUXiliari 

24 

ref gtr 

REFriGoraT R 

25 

r pht REsistanc SPace HaaT ng 

26 

r om c 

ROOM Air Conditioning 

27 

smidlm 

SMall InDustriaL H t r 

28 

wathtr 

WATor H aT R 


!(( )«> >t( 4t >*( tK 41 4i * >t( t >0 « % 41 *««« K( t * 


File 


load modelling 


NUMBER OF LOAD BUSES AT WHICH LOAD MODEL SIMULATION IS REQUIRED 
4 

Used with 


KV/P U) 


3E FREQUENCY 

Present frequency* 

Time step* ♦ 

30 0 

50 0 

0 05 

BUS NO 

REAL POWER 

BASEVOLTAGECof th bv 

8 

5 8 

132 00 

9 

11 2 

132 00 

11 

7 6 

132 00 

12 

22 8 

220 00 

14 

6 2 

1 00 

15 

8 2 

1 00 

17 

9 0 

1 00 

21 

17 5 

1 00 

30 

10 6 

1 00 


tability program 


*it.i|ii(.+)(i+i|ii((Xii|iH.i(i**)(c*)t)(i*Hi**#*****i(i*********i(i**'(<+**++*************>l'*+***** 


File 


load composition 


Total loads oi different models in MW 
Total number of models 20 AREA BUS NO 8 
MODEL NUMBER MODEL LOAD (in MW) 


1 

1700 

2 

0774 

3 

3190 

4 

0029 

6 

0029 

6 

0455 

7 

0661 

8 

1119 

g 

0029 

10 

0000 

11 

0081 

12 

4379 

13 

0461 

14 

0197 

IS 

0075 

16 

0044 

17 

0046 

18 

1 4947 

19 

3932 

20 

3492 

21 

1334 

22 

1331 
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23 


0625 




24 


14B0 




2B 


2378 




26 


2842 




27 


0196 




28 


8369 




LOSSES 



3832958 

MW 

T tal 

loads 

f diif r 

nt B 

dels 

in 

Total 

numb r 

of m dels 

28 

AREA 


MODEL 

NUMBER 

MODEL 

lOAD 

(in 1 

MW) 

1 


1 2800 




2 


1495 




3 


6160 




4 


0056 




S 


0056 




6 


4879 




7 


1277 




8 


2162 




9 


0056 




10 


0000 




11 


0157 




12 


8456 




13 


0890 




14 


0381 




15 


0146 




16 


0084 




17 


0090 




18 


9564 




19 


1 7248 




20 


6742 




21 


2576 




22 


2570 




23 


7000 




24 


2800 




26 


4592 




26 


2488 




27 


1241 




28 


8643 





LOSSES 7401676 HW 

and 80 on 


>|< i4( ito(i )K « itoti Hi 4< )ti ^ Hi ■t' If If >f * * <tc )ti 4i if 4i * ifi % )ti )*( If It* « i4< « Ht * ><( >t»|t >(i * * 


Fils 8 d 


STAriC AND DYNAMIC BEHAVIOUR PERCENTAGE 



NUMBER OF MODELS 




28 



MODEL NO 

STATIC 

DYNAMIC 

EXPONENTIAL/P OLYKDKIAL 




(stati ) 

1 

too 

0 

e 

2 

100 

0 


3 

100 

0 

e 

4 

100 

0 

6 

5 

100 

0 

6 

6 

loo 

0 

6 

7 

100 

0 

e 

8 

100 

0 

e 

9 

100 

0 

e 

10 

100 

0 


U 

100 

0 

a 

12 

100 

0 

a 

13 

100 

0 


14 

100 

0 

a 

15 

100 

0 

a 

16 

100 

0 

a 

17 

100 

0 


18 

100 

0 

a 


72 


19 

100 

0 

6 

20 

100 

0 


21 

100 

0 


22 

100 

0 


23 

100 

0 


24 

100 

0 

e 

2E 

100 

0 


26 

100 

0 

e 

27 

100 

0 


28 

100 

0 



ifiifif 4< 4 it"t< ><<>*< *l< >f >t< >!<>*■ >t> 4i If 4< 4< 4 4i * 4> 4< If * )t> % 1 ( 1 4: 4< >4i * 4> if ifififififif if>f >f>f )f’f4t>f>fif44<4i)f <f ifif 

Fil mp cl nr 

MODEL CHARACTERISTICS 
STATIC AND EXPONENTIAL 


NO 

cod 

por 

Pv 

Pf 

Ov 

1 

agriwp 

86 

1 4 

E 6 

1 2 

2 

arcfur 

72 

2 d 

1 0 

1 61 

3 

c ntac 

81 

20 

90 

2 2 

4 

cltdry 

09 

2 0 

0 0 

3 3 

E 

CltUBl 

OE 

08 

2 9 

1 6 

6 

cmceac 

76 

10 

1 0 

2 B 

7 

cmrmac 

76 

BO 

80 

2 6 

8 

colrtv 

77 

2 0 

0 0 

E 2 

9 

c oramhp 

34 

10 

1 0 

2 6 

10 

dshwsr 

99 

1 8 

0 0 

3 6 

11 

ellrya 

90 

1 8 

30 

2 2 

12 

flrlit 

90 

1 0 

1 0 

3 0 

13 

frezor 

84 

80 

60 

2 6 

14 

furfan 

73 

08 

2 9 

1 6 

IB 

hpcoac 

81 

20 

90 

2 6 

16 

hpcmnc 

81 

10 

1 0 

2 6 

17 

hpspht 

84 

20 

90 

2 6 

18 

inclit 

1 0 

1 64 

0 0 

0 0 

19 

Iridlm 

80 

06 

1 9 

60 

20 

ordf an 

87 

08 

2 9 

1 6 

21 

ordratr 

87 

08 

2 9 

1 6 

22 

ordpmp 

87 

08 

2 9 

1 6 

23 

pplaux 

80 

08 

2 9 

1 6 

24 

r fgtr 

84 

80 

60 

2 6 

26 

reapUt 

1 0 

2 0 

0 0 

0 0 

26 

rooraac 

76 

60 

60 

2 6 

27 

smidlm 

83 

10 

2 9 

60 

28 

waihtr 

1 0 

2 0 

0 0 

0 0 


Pf 


Nm 

POFnra 

P tun 

Pfnm 

Qvnm 

Qfnm 

4 

2 

1 0 

0 0 

9 0 

0 0 

0 0 

0 0 

1 

0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

2 

7 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

2 

6 

20 

1 0 

0 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

3 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

2 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

4 

6 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

3 

90 

1 0 

2 0 

0 0 

0 0 

0 0 

1 

4 

80 

1 0 

2 0 

0 0 

0 0 

0 0 

60 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

2 

8 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

4 

80 

1 0 

2 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

2 

7 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

3 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

3 

90 

1 0 

2 0 

0 0 

0 0 

0 0 

0 

0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

2 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

4 

80 

1 0 

2 0 

0 0 

0 0 

0 0 

0 

0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

2 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 

8 

1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 

0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 




File 


ouK plot 


To writs th« bua V 
Roqulremont of bus 

Raqulrsment of out 


P Q 


Number of buses 

Sorie.1 Number 
1 
2 
3 


if 

0 

Bua 


end total ayetem P nd 0 
result out put in a fileCYes 1 
1 

put files at all buses(Ye3 1 K 
0 

for limited buses is required 


Numbers 

1 

2 

3 


File Name 
0 1 


4 

11 

0 11 

6 

14 

0 14 

6 

17 

0 17 

7 

21 

0 21 

8 

27 

27 


in individual files 
Ho 0) 

0 ) 


I 



73 


9 30 30 

*♦♦*♦»**♦**** + ** ♦*♦♦»♦♦*♦♦*♦*** *%«*** ♦♦♦,*»**,.* 
File ut file namaa 

0 1 
0 2 

3 

4 

0 5 
0 6 

7 

and 3 on 

♦ ♦Hi*************************** *♦*♦** ♦***i(i«i|,********** *H.»»****» ***« 

File gen outag 

Requirement f grnorator outagaa or lose (Yes 1 N 0) 

0 

Number of gonaratoia 
3 

SL No Bus Itr Loss Factor 

1 4 6 0 0 

2 6 6 0 0 

3 6 6 0 0 

File load outage 

Requirement of load outages or loss (Yea 1 No 0) 

0 

Number of loads 
3 


SL N 

Bus 

Itr 

Loss Factor 

1 

14 

6 

0 0 

2 

27 

6 

0 0 

3 

30 

7 

0 0 


File itr limit 

LFP iteration limit 
20 

LFP+LHP iteration limit 
30 

Tolorenco limit of voltage for convergence 
0 001 

Requirement of varying tolerance limit of voltaga(Yes<"l H 0) 

0 

Niunber of variations 
2 

treation No Tolerance limit 

16 0 06 

26 0 01 

fffffeffifitigfififififififffiiffieififeeffeeeeeeeeifeeeeeeeeeeeeieee *♦♦♦♦♦♦♦**♦♦♦♦♦♦ 

File load mix change 

Requitement of load mix change (Yes 1 No 0) 

1 

Number of Iterations 
16 

m *iB«l **** ss a 


74 


ITR NO 1 

Numb r of bus s 2 

Bu and mod 1 number in increa ing order 
BUS NUMBER 8 
Numb r of m del 2 
Hod 1 6 Load change 0 0 

Model 28 Load change 0 0 

BUS NUMBER 11 
Numb r of mod la 2 
H d 1 6 Load change 0 0 

M d 1 28 L ad chang 0 0 

ITR NO 2 

Number of bus 2 

Bus and model numbers in incr aaing ord r 
BUS NUMBER 8 


Numb r 

of models 

2 


M d 1 

6 

Load change 0 

15 

Model 

28 

Load chang 

0 IB 

BUS NUMBER 11 



Number 

of models 

2 


Mod 1 

6 

Load hange 0 

IS 

Model 

28 

Load change 

0 IB 



a a 


and so 

n 




File react loss 


Requirom nt of mod lling of reactive loaa in diatribution(Yes 1 N 

1 

Number of buses 
9 


Bus no 
8 
g 

11 

12 

14 

IB 

17 

21 

30 


Reactive loss at 1 0 p u Voltage(MVAR) 
0 1 
0 35 

0 15 

1 2 
0 2 
0 3 
0 7 
0 3 
0 6 


0 ) 


If >(i >4< )<( )(oK Hi % iti ^ )4< )*( 4 ( If 4 * ifi 4 > )|i >4< >4* 41 « 


File load increment 


Requirement of load increment (Yes 1 No 0) 
1 

Total iterations 


IB 

Requirement 

Iteration 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


of pause(Ya8 1 No 0) 

0 

Factor of load increment 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 1 



75 


13 

14 

15 

T me rase 1 ad 

Numb X i buaoa 

Bus uum 
12 
21 


1 1 
1 1 
1 1 

at all bus (Yaa 1 No 0) 

1 

t which load is to bo mcraased 
2 


File d d 


MODEL CODE It NUMBER NUMBER OF RATINGS 

agriwp 1 2 


KW 

ISO 0 


NUMBER dynm modal numb r 
100 1 


1 0 1 0 0 05 
0 078 0 065 2 67 
04 0 01 0 01 


0 80 

0 044 0 049 

1 0 


KW 

300 0 


NUMBER dyn model number 
60 2 


1 0 1 0 0 05 
0 078 0 06S 2 67 
04 0 01 0 01 


0 80 

0 044 0 049 

1 0 




MODEL CODE ft NUMBER 
centao 3 


NUMBER OF RATINGS 
2 


«»**%* iti If **** Hi * 1*1 *<)"•< If iH * 't'*** ^ If **** 'tot< 4> **** 


KW 
2 5 


NUMBER dyn model numb r 
3000 3 


1 0 1 0 0 05 
0 078 0 066 2 67 
04 0 01 0 01 


0 80 

0 044 0 049 

1 0 


KW 
5 0 


NUMBER dyn modal number 

1100 4 


1 0 1 0 0 05 
0 078 0 06B 2 67 
04 0 01 0 01 


0 80 

0 044 0 049 

1 0 




File itr outl 

In this the voltages real and reactive powers at all buses 
and total real and reactive power in the system and los es 
are written ior each load increment iteration 


- 


ITR No 10 

- 

Base HVA 

•> 100 0 

Pgen 

2 9461 

Pload 


2 79046 

LOSaQB 

» 1556332 

Bus N 


V 

Pinj 

qioad 



Cp u) 

(p U) 

(p u) 

1 

1 

0600 

1 5498 

4258 

2 

1 

0333 

6981 

1200 

3 


9859 

6981 

1000 

4 

1 

0259 

0000 

2400 


76 


5 

1 0100 

0000 

1528 

6 

1 0008 

0000 

2400 

7 

9773 

0000 

0000 

8 

9478 

1168 

0329 

9 

9672 

2325 

0981 

10 

1 0357 

0000 

3790 

11 

9767 

1576 

0487 

12 

9704 

4751 

2145 

13 

9748 

0000 

0000 

14 

9361 

1238 

0370 

15 

9307 

1630 

0507 

16 

9446 

0698 

- 0359 

17 

9372 

1840 

0813 

18 

9090 

0638 

0180 

19 

9032 

1895 

0678 

20 

9125 

0439 

0140 

21 

9290 

3559 

1555 

22 

9314 

0000 

0000 

23 

9213 

0638 

0319 

24 

9276 

1735 

0479 

25 

9830 

0000 

0000 

26 

94B6 

0698 

0459 

27 

9913 

0479 

0239 

28 

9735 

0000 

0000 

29 

9911 

0479 

0180 

30 

9636 

2119 

0561 


ITR 

No 11 

- 
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File itr out2 

In this the voltage real and reactive pow ra at the modelled buses 
IS written for each LFP+LMP run 

ITN BUS NO VOLTAGE P q 

CKV) (MW) (WAR) 


BASE LOAD CONDITION ♦♦*♦♦♦♦♦♦♦»* 



(LFP 

+ LHP) ITR 

1 - - 


1 

8 

1 0347 

6 0187 

1 9742 

1 

9 

1 0493 

11 6307 

6 6877 

1 

11 

1 0247 

7 8117 

2 6966 

1 

12 

1 0085 

22 9410 

11 4775 

1 

14 

1 0351 

6 4268 

2 2304 

1 

15 

1 0319 

8 4680 

3 0893 

1 

17 

1 0298 

9 1822 

4 8043 

1 

21 

1 0252 

17 8018 

9 2423 

1 

30 

1 0614 

11 0657 

3 3747 


(LFP 

+ LKP) ITR 

2 — - 


2 

8 

1 0385 

6 0304 

1 9816 

2 

9 

1 0501 

11 6378 

5 6966 

2 

11 

1 0240 

7 8065 

2 6926 

2 

12 

1 0078 

22 9294 

11 4640 

2 

14 

1 0346 

6 4237 

2 2286 

2 

15 

1 0320 

8 4686 

3 0897 

2 

17 

1 0321 

9 1963 

4 8235 

2 

21 

1 0276 

17 8306 

9 2806 

2 

30 

1 0413 

10 9330 

3 3068 


(LFP 

+ IMP) ITR 

3 __ 


3 

8 

1 0366 

6 0302 

1 9814 

3 

9 

1 0501 

11 6378 

6 6966 

3 

11 

1 0241 

7 8068 

2 6927 

3 

12 

1 0079 

22 9299 

11 4646 

3 

14 

1 0346 

6 4238 

2 2286 

3 

15 

1 0320 

8 4686 

3 0897 

3 

17 

1 0320 

9 1960 

4 8232 


3 

21 

1 0276 

17 8301 

9 2799 

3 

30 

1 0421 

10 9439 

3 3114 

♦+**♦*♦*♦♦ LOAD 

INCHEHENT 

ITR 1 


(LFP + 

LHP) ITR 

1 


1 

8 

1 0316 

6 6417 

2 1879 

1 

9 

1 0468 

12 8811 

6 2293 

1 

11 

1 0212 

8 6186 

2 9482 

1 

12 

1 OOB2 

28 3396 

12 8836 

1 

14 

1 0297 

7 0789 

2 4292 

1 

IB 

1 0266 

9 3261. 

3 3646 

1 

17 

1 0267 

10 1448 

8 2864 

1 

21 

1 0217 

19 6622 

10 1042 

1 

30 

1 0369 

12 0422 

3 6038 


(LFP + 

LHP) ITR 

2 


2 

8 

1 0318 

6 6432 

2 1889 

2 

9 

1 0469 

12 8823 

6 2306 

2 

11 

1 0213 

8 6189 

2 9484 

2 

12 

1 0082 

28 3402 

12 8843 

2 

14 

1 0298 

7 0772 

2 4300 

2 

IS 

1 0268 

9 3279 

3 3688 

2 

17 

1 0269 

10 1463 

8 2888 

2 

21 

1 0219 

19 6686 

10 1087 

2 

30 

1 0373 

12 0474 

3 6068 


)ti id «««»«>(( load increment itr 






